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ABSTRACT. The thermal conductivity of Cu, Ag, Au, Mg, Cd, Al, In and one poly- 
crystalline and two single crystal specimens of Zn have been measured in the range 2—40° K. 
Each specimen shows a maximum in the conductivity in this range. The values of the 
coefficients « and f# in the equation for the thermal resistance R=aT7*-+-£/T have been 
calculated for each specimen. 


$1. INTRODUCTION 
ERY few data are available regarding the heat conductivity of metals at 
\ / low temperatures. Most of the experimental work has been done on 
superconductors and on the change in heat conductivity during the 
superconducting—normal transition (e.g. Mendelssohn and Olsen 1950, Hulm 
1950). Except for one or two isolated metals and alloys we have little information 
on the non-superconductors. In order to be able to test current theories and 
also to provide’ information for their improvement, it was decided to study 
systematically the thermal conductivities of the metallic elements. Since even 
minute amounts of impurity can cause an appreciable change in the heat 
conductivity, only specimens which were in a high state of purity were used. 
In the main these were Johnson Matthey spectrographically standardized 
substances. 

The specimens were in rod form about 5cm long and, as a rule, between 
1 and 2mm in diameter. A small electric heater was wound around one end and 
the other end was connected to the expansion chamber of a Simon helium 
liquefier. Helium gas thermometers were attached to the specimen by copper 
contacts 3cm apart. The pressure difference in the thermometers corresponding 
to their temperature difference was read on a differential oil manometer. As a 
rule measurements were taken between 2 and 40°k. ‘The relative accuracy for 
the results for one specimen was about 2%, except for very good conductors, 
where the temperature difference obtainable was very small. ‘The absolute 
accuracy was usually limited by the determination of the form factor of the 
specimen and was of the order of 2 to 3%. A more detailed description of the 

experimental arrangement and procedure will be published later. 
This first paper gives the results so far obtained for the metals of groups 1, 
2 and 3. A further paper will deal with those transition elements that have been 
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investigated so far. Since these communications are both in the nature of 
preliminary reports, no detailed analysis of the results is presented. 

The metals used are shown in table 1. The annealing was done in evacuated 
quartz tubes for several hours at two-thirds the melting temperature. 


Table 1 

Origin Purity % Heat treatment and crystal structure 
Copper 1 JM 4234 99-999 Annealed polycrystalline 
Silver 1 JM 1722 99-99 Polycrystalline wire 
Gold 1 JM1916a 99-999 Polycrystalline wire 
Magnesium ‘1 JM 1703 99-95 Polycrystalline 
Zinc 1 HS 8392 99-9995 Polycrystalline 
Zinc 2 SC: 99-997 Single crystal; hex. axis 80° to rod axis 
Zinc 3 ESIC: 99-997 Single crystal; hex. axis 13° to rod axis 
‘Cadmium 1 HS 10547 99-9999 Cast in glass 
Aluminium 1 JM 4899 99-994 Annealed polycrystalline 
Indium 1 JM 4398 99-993 Polycrystalline wire 


JM=Johnson Matthey. 
HS=Hilger H.S. brand. 
1.S.C.=Imperial Smelting Corp., Avonmouth. 


§2. RESULTS 

The results are shown in figs. 14. In all cases the conductivity has a 
maximum in the region investigated, and the value of the conductivity at this 
maximum is of the order of 10 watts/cm/deg, except in the case of copper. 
Electrical resistance measurements on all the specimens have not yet been made, 
but it seems unlikely that there is any appreciable lattice heat conduction. Hence 
we assume that the thermal resistance R can be written in the form 
R=aT?+/T, where « is a term depending on the scattering of the electrons 
by the lattice vibrations and 6 is an impurity scattering term. In each case 
graphs of RT against T* have been plotted and the coefficients « and 8 have been 
determined for the low-temperature part of the curves, since in no case are they 
linear over their whole extent. ‘These values are shown in table 2. 


Table 2 
Specimen Cull "Agd “Au Me te7n And) sone Cda Ala maetriat 
ax 105 392, 9:07 18 10:6 21 34 31 140 2-2 190 
B O:35)" 1:60) Ai 1>25y s0548 2 0-7eeeen 0 GnenO 2 5am SO 


Copper, Silver and Gold 

The values of the thermal conductivity are shown in fig. 1. The copper 
specimen was rather too thick to measure the extremely high conductivity, and 
this accounts for the scatter in the points. Berman and MacDonald (1952) have 
obtained a thinner specimen from the same batch of copper, and _ their 
conductivity curve has a much higher and sharper maximum at approximately 
the same temperature as ours. The curves agree at the higher temperatures. 
The higher maximum is presumably due to the difference in the annealing and 
preparation of the specimens. The conductivity of the gold and silver is much 
less than that of the copper. The broad maximum of the silver seems to indicate 
a large impurity scattering term, and it is hoped to repeat the determination with 
a purer sample of the metal. 
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Zinc, Cadmium and Magnesium 
The results for the three specimens of zinc are shown in fig. 2. Zn1 is a 
polycrystalline specimen. Zn2 and Zn3 are single crystals from the same batch 
of metal, but differing from Zn 1. They have their hexagonal axes at 80° and 13° 
to the rod axis respectively. These directions were the nearest that could be 
obtained easily to axes parallel and perpendicular to the rod. It will be seen 
that the curves for Zn2 and Zn3 coincide at high and low temperatures but 
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cadmium and magnesium, aluminium and indium. 


that Zn3 has a maximum about 10° greater than that of Zn2. ‘This probably 
has some connection with the fact that the atomic spacing in zinc is different 
perpendicular and parallel to the hexagonal planes, but experiments on more 
single crystals need to be done in order to clarify this point. 

Figure 3 gives the conductivity curves for magnesium and cadmium. ‘The 
cadmium specimen is the one which showed a large increase in thermal resistance 
in a magnetic field, as has been reported previously (Mendelssohn and 
Rosenberg 1951). ; 
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In metals such as gold and magnesium, where a minimum often occurs in 
the electrical resistivity, no such corresponding effect was found in the thermal 
resistance. The electrical minimum is usually very shallow, and in the thermal case 
this would be reflected as only a very slight change in the slope of a curve, and 
not, as in the electrical case, by a change from a constant value. It is not 
surprising, therefore, that no analogous results to the electrical minimum have 
been found in these experiments. 


Aluminium and Indium 

The conductivity curves are shown in fig. 4. The indium showed an 
appreciable increase in thermal resistance in a magnetic field, and the values in 
the normal state below the superconducting transition temperature have been 
obtained by extrapolating to zero magnetic field. It is interesting to compare 
these results with those of Hulm (1950), whose indium specimen was less pure 
than ours and showed no maximum. The value of « obtained is the same as his, 
which should be the case if the simple equation for R is correct. On the other 
hand, the results for aluminium do not agree with those obtained by Webber, 
Andrews and Spohr (1951) whose value of « is 2°7x10-°. Their specimen, 
however, may be slightly purer than ours. 

The effect of a magnetic field on the thermal conductivity of these specimens 
has been determined and will be published separately. 
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ABSTRACT. ‘The thermal conductivity in the range 2-40°K has been measured for plate 
Mn, Fe, Ni, Zr, Cb, Mo, Rh, Pd, Ta, W, Ir, Pt, U, and for a Pb single crystal. Asa rule, 
the metals of group 8 have a much higher conductivity than have the metals of the earlier 
groups. ‘The values of « and f in the equation for the thermal resistance R=aT?+ B/T have 
been calculated for each specimen. The experiments on Cb indicate that the normal state 
was not reached in a field of 3300 gauss at temperatures below 5°K. 
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§1. INTRODUCTION 


N a previous paper (Mendelssohn and Rosenberg 1952) we have dealt with the 
| thermal conductivity of some metals of groups 1, 2 and 3. This paper 

presents the results of the rest of our work so far completed. In addition 
to some transition elements, it also gives the results of experiments on a lead 
single crystal and on uranium. The specimens used are described in table 1. 

Experimental procedure is similar to that outlined in the paper referred 
to above. 

§2. RESULTS 

The graphs of thermal conductivity are shown in figs. 1-7. The most 

outstanding result is that the conductivity of the metals at the beginning of each 


transition group is very small compared with those at the end of the group, 
ie. the group 8 elements. For those elements with low conductivity no 


Table 1 

Specimen Origin Purity (%) Notes 
Titanium 1 A.E.I. 99-9 
Titanium 2 ALE, 99-99 Annealed 
Manganese 1 JM 2472 99-99 Annealed 
Iron 1 JM 4975 99-99 Annealed 
Nickel 1 JM 4884 99-997 Annealed 
Zirconium 1 JM 3062 ~98 Annealed 
Columbium 1 JM 4526 99-99 
Molybdenum 1 JM 2331 99-95 
Rhodium 1 JM 2357 99-995 Wire 
Palladium 1 JM 2134 99-995 Annealed wire 
Tantalum 1 JM 3804 99-98 
Tungsten 1 JM 2260 99-99 Annealed 
Iridium 1 JM 3441 99-995 Annealed wire 
Platinum 1 JM 2157b 99-999 Annealed wire 
Lead 1 Tadenac 99-998 Single crystal 
Uranium 1 ACHR. 


JM=Johnson Matthey, A.E.I.=Associated Electrical Industries Research Laboratories. 


Table 2 
Specimen «105 B Specimen ax 10° B 
Anya — 290 Rh 1 22 1-4 
Abe Ps — 170 Pdi 64 11 
Mn 1 “= 1200 wart = 27 
Fe 1 18 9°5 Wil 10:2 5:9 
Ni 1 22 4-4 tea 3-6 0:77 
Zr 1 130 76 ipre al 43 0:40 
Cpe — — Pb 1 325 0:10 
Mo 1 75 6:7 et 750 95 


maximum was found in the experimental range covered. In some cases the 
curves were linear within the limits of experimental error, and values of the 
coefficient « in the equation for the thermal resistance R=a7T?+8/T could not 
be evaluated. In all other cases graphs were drawn of RT against T? and the 
values of « and £ for low temperatures were calculated. ‘These values are shown 
in table 2. 
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Titanium, Manganese, Iron and Nickel | 

The results for these metals are shown in figs. 1 and 2. The manganese — 
and titanium both have very small conductivities of the form dictated by impurity 
scattering. On the other hand, the low conductivity of manganese may have 
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Fig. 1. The thermal conductivity of uranium, titanium and manganese. 
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Fig. 2. The thermal conductivity of nickel, iron, molybdenum, palladium and tungsten. 


some connection with its extremely complicated crystal structure. Ti2 has 
been annealed and may be slightly purer than Til. Nickel and iron show the 
large increase in the conductivity for the group 8 elements. 


| 
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Fig. 3. The thermal conductivity of zirconium and tantalum. 
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Fig. 4. The thermal conductivity of columbium. 
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Zirconium, Columbium, Molybdenum, Rhodium and Palladium 

Zirconium (fig. 3) and columbium (fig. 4) again have low conductivities 
followed by the higher values of molybdenum (fig. 2), rhodium (fig. 5) and 
palladium (fig. 2). 

Columbium being a superconductor, a magnetic field has to be applied in 
order to measure the heat conductivity in the normal state. The curves show 
that this state is not reached at low temperatures even in a field of 3300 gauss, 
which is much higher than the published values of the transition field. Assuming 
that the conductivity curve in the normal state is a straight line through the 
origin, then the 2300 gauss curve breaks away from this about 6°K and the 
3 300 gauss curve at a slightly lower temperature. The 2300 gauss curve crosses 
the superconducting state curve at about 3°K. One explanation for this is that 
the lattice conductivity is decreased because of the extra scattering by the normal 
electrons produced by the field, and this decrease is not fully compensated by 
the increase in the electron heat transport. Electrical measurements on this 
specimen, however, indicate that almost all the heat transport is by electrons, 
and it may be that it is the electron heat transport itself that is decreased by the 
magnetic field. A possible mechanism is by the formation of normal and 
superconducting laminae, the boundaries of which may tend to scatter the 
electrons. However, further experiments, when larger fields are available, will 
be necessary to investigate the superconducting transition of columbium properly. 

Zirconium usually has a large impurity content of hafnium, and this might 
- account for the low conductivity found. On the other hand, the very great 
similarity in structure and properties of the two metals, as is shown by the great 
difficulty presented in their separation, may not cause the conductivity to be 
reduced as much as one would otherwise suppose, i.e. the hafnium might not 
act as a true impurity. Molybdenum in group 6 has a much higher conductivity 
whose value is nearly the same as that of the next group 6 metal, tungsten. 
Rhodium and palladium both have conductivities of the order of a few 
watts/cm/deg like the other group 8 metals investigated. The palladium shows 
a flat maximum at about 24°. 


Tantalum, Tungsten, Iridium and Platinum 

The curve for the normal state of tantalum is shown in fig. 3 and the curves 
below the superconducting transition temperature are given in fig. 6. In this range 
the normal curve does not seem to be strictly linear. 

Tungsten (fig. 2), referred to in the preceding section, has a much higher 
conductivity, with a maximum at about 28°. The value for « calculated by 
Hulm (1950) from the results of Bremmer and de Haas (1936) is 0-88 x 10-4, 
and it compares fairly well with our value of 1-02 x 10-4. 

Iridium and platinum (fig. 5) have high conductivities, that of iridium being 
the highest of the transition metals so far measured. Both metals exhibit maxima 
in the range covered. 

Lead 


The presentation to us by Mr. S. Weintroub of a lead single crystal gave us 
an opportunity of measuring its thermal conductivity in the normal and 
superconducting states. ‘he specimen was extremely pure, as is shown by 
its very high and sharp maximum in the normal state (fig. 5). Although readings 
were taken on the low-temperature side of this maximum in the region of 1-8°x, 
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the conductivity was changing so rapidly that it was extremely difficult to get any 
reliable experimental points in this range, and hence only one point is shown on 
the graph for this side of the curve. The heat conductivity in the superconducting 
state is shown on a larger scale in fig. 7. This confirms the sudden breaking 
away of the superconducting from the normal curve and also the shallow minimum 
_. in the superconducting curve at about 5°K, as has been found by other workers 
(de Haas and Rademakers 1940). The rise in conductivity on the low-temperature 


20 5 
4; 
S 
» 
n=] 
es: 
a FI 
BS 3 
= 0 
= 
= = 
> _ 
c = 
S Ss 
ue 
0 10 20 30 40 L 
Temperature (°k) 0 2 Aerys. tu § 8 
Fig. 5. The thermal conductivity of iridium, Temperature (°K) 
rhodium, platinum and a single crystal of lead Fig. 7 The thermal conductivity of lead in 
in the normal state. the superconductive state. 


Normal 


Superconducting 


Conductivity (W/cm/deg,) 


ar 1 4 1 
2 25 30 35 40 45 
Temperature (°K) 


Fig. 6. The thermal conductivity of tantalum in the normal and superconductive states. 


side of this minimum has a connection, presumably, with the extremely rapid 
increase in the normal conductivity in this temperature range. 

Magneto-resistance effects, which will be reported separately, have been 
found in this specimen. 

Uranium 

The curve for uranium is shown in fig. 1. Although it has a low conductivity 
it is not linear, whereas all the other metals investigated which have a conductivity 
of this order are linear in the low-temperature range. Electrical resistance 
measurements indicate that, within the experimental error, the whole of the heat 
transport in the uranium is by electrons. 
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As has been mentioned previously, this paper is only in the nature of a 
preliminary report, especially as regards the analysis of the results. ‘The graphs 
of RT against T? are by no means linear even at temperatures where the simple 
formula for R should hold, i:e. below one-tenth of the Debye @ temperature 
(Makinson 1938), and it is clear that the details of the theory need revising in order 
that our results may be explained satisfactorily. Extension of the experimental 
work to higher temperatures is in progress. 
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ABSTRACT. The vibrational frequencies of an orthorhombic ionic lattice, of the 
thallium fluoride type, have been worked out in detail for wavelengths long compared with 
the crystal spacing though short in relation to the crystal size. It is shown that these 
frequencies are not unique but depend in general on the direction of the wave in the crystal— 
in contrast to what happens in cubic crystals. The effect of this behaviour on the vibra- 
tional spectrum is discussed briefly. 


§1. INTRODUCTION 

T has been shown by Lyddane and Herzfeld (1938) and by Kellermann (1940) 
| that the normal modes of an alkali halide (of the rocksalt type) of long 

wavelength and high frequency consist of two transverse vibrations of 
equal frequency and a longitudinal vibration of much higher frequency. The 
term Jong wavelength here means long compared with the crystal spacing, though 
short compared with the size of the crystal. Further, the frequencies are 
calculated on the assumption that Coulomb forces (as well as repulsive forces) 
would suffice; this excludes cases where the crystal radiates electromagnetic 
waves—at any rate pending a special consideration—but, as far as those 
frequencies are concerned which play an important part in the vibrational 
spectrum, this restriction is not significant. 
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The existence of the two types of limiting frequencies is linked on the one 
hand with the existence of a long high-frequency tail to the vibrational spectrum 
(composed of mainly longitudinal vibrations), on the other hand with a large 
peak in the vibrational spectrum (associated with mainly transverse vibrations), 
as can be seen from the work of Kellermann and of Iona (1941). 

The large difference in the frequencies of the transverse and the longitudinal 
waves is a characteristic ionic effect. This effect is considered in some detail 
in the present paper, particularly for crystals which depart from cubic 
symmetry. Here a new aspect is found. If a lattice is considered containing 
two oppositely charged particles in the unit cell, three non-zero frequencies are 
always obtained in the limit of long waves; if the lattice has, for example, different 
properties along its axes (e.g. different spacings), some of the solutions of the 
equations for the frequency would be expected to be different if cases are 
compared where the wave vector is chosen along different axes. For instance, 
a purely longitudinal wave, of long wavelength, should have different frequencies 
according to the axis along which the wave vector is directed. 

This suggests that the frequencies are not unique in the long-wave limit 
except in special cases such as that of the cubic crystal. As variability in the 
frequency would have an effect on the vibrational spectrum, which is unexpected, 
a detailed study was made of a special case. 


§2. THE ORTHORHOMBIC CRYSTAL 

The case of an orthorhombic crystal seemed particularly suitable for this 
investigation. ‘The particular model chosen was based on thallium fluoride, 
which contains four molecules per orthorhombic cell. ‘This cell is similar to 
that of rocksalt and can be produced from it by appropriate dilatations along 
the orthogonal axes. Furthermore, a smaller cell containing two particles can be 
chosen, as suggested by the rhombohedral cell in rocksalt, and this cell will 
be used in the calculations below. The definition of the lattice vector also 
follows the same pattern as in rocksalt, and the mathematical descriptions and 
formulae can always be reduced, for the purposes of checking results, to that of 
the rocksalt case. 

The lattice vector of the cell r can be written as 


pals = La, = AEN Sie Las, eoveee (1) 
where a, =61,+cl,; a,=Cig+ai;; ag=ci, big. ss a ass (2) 


Here i,, in, i, are unit vectors along the x, y, z directions respectively, these 
being parallel to the orthorhombic axes. The smallest distance between similar 
particles along these axes are respectively 2a, 2b, 2c, where in thallium fluoride 


Zao BAT 2b=5:50'4,2Ze=6-08 Ke (3) 
The particle vectors r, in the basic cell and the charges «, are given by 
r,=0, r,=43(a,+a,+4s) 
Cpe | Cl tes ot |G ee Me 1 endo nsis (4) 


where e is the electronic charge. 

The vibrational frequencies of long wavelength can be obtained as a special 
case in the general calculation of lattice frequencies. ‘The method of doing this 
is well known (Born 1923). Taking a particular wave of wave vector s, 
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wavelength A, and angular frequency w, the displacement u,' of a particle k in 
the /th cell, where /=(/,/,/;), is given by 


u,/=U, exp (—zwt) exp {2mi(r,'.S/A)}. eee ee (5) 


The usual boundary conditions (of the periodic type) will determine the possible 
values of s/A. 


The equations of motion are 


wm,U yet > » os ] Uy, =9, o Tonsnainelia (6) 
key L*Y 


where ese =a (bie, exp 2at( Ky S/A) on eet eee (7) 


Here ri,,=ri, —r,° is the vectorial distance between particles of type k’ and 
type k, the first being in the /th cell, the second in the reference cell /=0. The 
quantity ¢/,;, is the potential energy of this pair of particles, and the suffixes x, y 
denote differentiation with respect to these directions. The masses of the 


particles are given by m,. Further, in the calculation of ee the terms 


(P17) xy are determined by the condition 
2 (Peleg =ULt oe) ae Or ee (8) 


As there are two particles per cell, the angular frequencies can be determined 
from a sixth-order determinant (cf. eqn. (6)), once the lattice sums (7) have 
been evaluated. For the special case of long wavelength, three of the solutions 
will be zero; hence the determinant is reducible to the third order, or even 
further if waves are chosen in special directions. 


§3. CALCULATION OF ELECTROSTATIC LATTICE SUMS 


In the calculation of the lattice sums, the electrostatic part needs special 
methods; equivalent methods are known, one due to Ewald (1938), the other 
using Epstein zeta functions (Epstein 1903). The latter method has been used 
by Broch (1937) in treating the linear ionic lattice, the former by Kellermann 
(1940) in connection with the vibrational spectrum of rocksalt. As the 
Epstein functions seem to lead to a more direct calculation and deserve to be 
better known, they have been preferred here. 

The function which we shall use is defined as 


g , exp (271 X Lh,) 
Z => ———__F* , BS cata Ge, 
n |) == Tare ©) 
where (l+g) = > Cae +2,)(1,+2,) 


ae 
and yw, ¥=1, 2,3. The C,,, are constants, the J, are integers ranging from — o 
to + 0. ‘The h, range from 0) to 1, and sis an integer. The denominator in (9) 
is in the form of a power of the lattice distance (|r4,,|)* where C,,=C,, and 
the g, are numerically equal to the coordinates of the particles in the cell. 
If, further, we put 4,=(a,.s/A), then the zeta functions, suitably differentiated 
with respect to the h,, will be in the form required for the evaluation of the 
electrostatic part of the lattice sums (7). To achieve this evaluation numerically 
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we need the zeta function in a quickly convergent form, which has been given 


by Epstein (1903): 


-on() 


exp {—271 (asa 
+ RU | _ ds 20-908 exp (—n2@(h-+m)—2ni ¥ mg, 


(s) = 6 dz 2-98 ¥ exp {—nayi(g +) +2ni XI} 


Here I’ refers to the gamma function and A is the determinant formed by the 
C,,,. The quantity ¥(g+/) has been given above, and ®(h +m) is defined by the 
following expression : 


Diem) =" Coie an (he m,), (by v= 1, 263) nel oys (11) 
“uy 
where the m,, are integers ranging from — 0 to +o, and 
mee 21 TON 
Cy = ADO ee es (12) 


‘The transformation (12) is a reciprocal one, and it can be shown that O(m) is the 
square of the distance in the reciprocal lattice if »s(/) is the square of the lattice 
vector. ‘I'he summation over the integers /,, m, extends also from — oo to + ©. 
In our calculation p=3. 

The integrals in (10) could be used in a slightly different form by the 
introduction of the Ewald separating factor «9, see Born (1923), instead of taking 
€)=1, but it has not been found necessary to introduce this modification here. 


§4. DETAILED EVALUATION OF THE ELECTROSTATIC 
LATTICE SUMS 


: 12m : : 
We consider first the evaluation of Fa . The lattice distance |7,,' | is 


given by 
(ra!) =20°[$1,2(a +B) + HL +B) + $y2(1 +a) + ly +a ly 
+ BL, —-L(«a+f)—1(1+f8)—-L(1+«)+4+4a+36]. «0... (13) 
where CD ye Cd eee © ©) Ne er stcee. (14) 


It is found convenient to use the Epstein formulae in a dimensionless form, 
and hence we choose ¥(/+g) equal to the quantity in the bracket of (13). Hence 
the coefficients C,,,, g, are given by the following scheme: 


Cy =3(%+B8); Cr=Cy =38; 21= 


Coo=3(1+8); Cig=Cy=30; = 
Csg=3(1 +a); Cos=Cyo=3; 83= 33 


Die 


| 
Die ble oe 
~ - 
——— 

. 

. 

. 

. 

. 

. 

— 

— 

on 

N 


and A=4(«8). We can further determine the coefficients C’,, using (15). This 
leads to the following set of relations: 


, / 1 a aE ! / oe 
C1 = Cog = C53 = : a C2 = Cy = a 
ae (16) 
J 


, , =e ae , r Che a eat 
Cp= Cin = ES Ce geet ae ace 
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Finally we can write down ®(h + m) in terms of the C’,,,: t 
O(h +m) = Sr th 24+ My? + hz? + m,? + m2 + mg") 
re <s LE Chal wae ane | 
basses Ponte “i cham icley (9 ety is 2 | 
ae bare ae Bp (Aphz + mgm + hym3 + hgmy) 
fh ies ere: OF (hsm, + liga a) (17) 


With these preliminaries over, it is now necessary to consider how the 
derivatives of the potential occurring in the lattice sums can be formulated in 


: a 
terms of the zeta function. In fe we need the term 412(|742'|). Now 


19 Ole, = ep) Se Ae (18 
io [712° P? ee Ce vox Sa lis CR 
where 
Aye =[— («+ B)L? + (2 — By? + (2 — alg? + 441, — 2acl;ls — 28h, 
+2(a+ B)l; -—2(2—8)l,-—2(2—a)l,+(2—-a—f)]}.  ...... (19) 
Hence (as s=5) 
Wass e : F F 
ea = — 9573 EXP {mi(h, +hs+hs)} | - Z ; + 5 +Q-a-p)F, | 3 
eee (20) 
where 
OZ OZ eZ OZ CA OZ 
F,= <(2-—«) -—. +(2— 8) — — — +45-—— —2«—— — 28 —— 
\( “ape *C—Plaps —(+Blapa + tap ap — 2p ap. 28 age 


| QZ AZ aZ 
Ea are ee oa 
F,=Z. : : 


Inserting the values for the derivatives of Z, using (10), we obtain 
12 eiles A a 
ae She are PE dz 20° 5 HLE(L) exp (2ni% (l,—3)4,) 


xexp {—m2¥(I-1)} — a [° _ 4% G,,(mh) 


xexp {—73O(h+m)+71 & m,3| To, ee tel, Prete (21) 
where H7; has been defined above (eqn. (19)) and 
G,,(m, h) = a 2+ hy? + hy? + m,2 + m2 + m2) 
— re (Ayhy + mymy + hm, + hym,) 
— eee (Ayh; + mms + hms + hm) 


2(2a8 + 0% + 
Si <oreete) (Ash; + MgMs + hams + hymyg). Pits: (22) 
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The expression for G,,,, is independent of k, k’, and so these indices are omitted. 
We shall be concerned only with those cases ‘where the h, tend to zero. This 
is not necessarily the same as putting /,,=0, but is, in fact, the same except for 
the one term in G,, where m,=0. ve term is handled separately, i.e. 


ql 
ie wee Ne ot 2 emi 1(5/2) A vz lim | dz G',,(h) exp {—72z®(h)}, 


where ®(/), G,,(2) are obtained by putting m,=0 in (17) and (22) respectively. 
In the first term on the right-hand side all #,=0; in dealing with this term 
numerically it is convenient to transform the summations so as to make use 
of the Boe cell. Then 


ql? 7 00 
22 | 12 
is Biker -< 2527 5/2) 21 (5/2) | ' dz 23? a Gil P) 
iI ote) 
x exp {—4na(p.2-+ aps" +B) ~ am], 4 = Bax(m) 


x exp { — 3m (ont + ee a “) + at(m, + mM, + m), .. (24) 


where in the first summation the integers p=(P,,~2,p~3) range over the 
(orthorhombic) face-centred lattice of the particles k=2, the origin being taken 
at one of the particles k=1. The second summation m=(m,, mz, ms) is taken 
over the orthorhombic body-centred reciprocal lattice. The expressions for 
G2, B,, are given below. 

The second term in (23) can be evaluated directly, but the result is of interest 
only when special values are chosen for the h,. For instance, if h; =h,=h and 
hz=0, Gez,= —4h?/B and O(h)=2h?/8. ‘The term involving the integral in (23) 
then becomes 

dvet yrtleee 12 Aq Cte 
Zea V(5{2) AM? 3 2a? («p)t2 
since [(5/2) =321#/4 and A=4af. 

This term is typical for Coulomb forces, and for «=1, B=1 reduces to the 
value which appears in the calculation of the limiting frequencies of rocksalt. 

The other lattice sums of the above type can be obtained using the following 
list : 

1 a 
CH= (2p,2—aps?~Bpat)s Bea= (2m2-+ =m? — 3m?) 
Z 1 
CH= (Zapp? — Bp,”—p,"); Byy= € mM,” — B ms” — m:) 
j} 1 
C2 =(2Bp3"—p,*—ap2"); B= (; ms” — mm,” — e m2) 
(28 — a —«B) 2(28 + % — af) 
G,,(4) = rs pana (1? + hy? + hs”) — SS Sah hyhy 
2(2B+a+a 2(a —28 —«B) 
oe Bae) eR a hah, 


as 
bo 
ron’ 
— 


ne aera i 


Aral: B—2«) 2(B —2«— oy kes 
ieee hyhg+ ae ae | 
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The next set of lattice sums which need to be considered are of the type 


ss E 
el , xy. Here it can be shown that heal =( because of the symmetry 


h=0 
properties of the orthorhombic lattice; it is symmetrical on reflection in the 
vy, yz, zx planes. Using the same methods as in calculating (18), 


EK ac) 
a = = ayyA lim | Gay(h) exp (720 (H)}, oeee (27) 
XY Nn->0 h>0/ 1 
3@n are 1. eatieye 
— @ 2°1(5/2) Atle = we (ap)? Cae en YC fh 
The other quantities are given below together with their forms in the other cases: 


where 


Agy =2385 Gay () = = (ig? — fy? hy + Zhyha) 
Cpt Gh) = 3h — he 3? — hy? + 2hyhs) Ps cane (29) 
aye = (#8); Gy) = - 5 (4:2 — ig? — igh + Dlg). 


So.) LAE COMPEE TLE, DAL PICE SUM 
The contribution from the repulsive forces is considered here in its most 
elementary form in order to show the effect we wish to demonstrate in its most 
striking form. We assume that these are nearest neighbour repulsive forces, 
and further that these are independent of the direction of the bond. The more 
natural assumption of different force constants for different directions will 
merely introduce extra, and unknown, parameters. From the above assumptions 


it follows that 
12]® 12]® 1278 2 
Bae ee le Ps 
h=0 h=0 SS |p=0 3 2a? 


Ee Sale =[ seh 7° | ay (30) 


where f is a numerical Ais The evo tre form chosen for D is convenient 
in that this form occurs naturally in the electrostatic terms. The value of f is 
chosen so as to bring the frequencies considered into the same region as that 
of rocksalt (f=1-88, cf. §6). There is also in (30) no difference between h->0 
andh=0. The sorniee lattice sums can therefore be written as 


(ane ergo 
[: alt oe apd eee i ae (31) 


11 
The calculation of the lattice sums of the type | ] is carried out in the same 
Ng 


manner as the type we have been discussing, except that there is an extra term 
in the summation, involving (¢/,),,,-. This has to be determined using (8). 
It can then be shown, tor instance, that 


[Ler es “Leon” Letfnw Les 
XX h>0 5.064 h>0 ieee Bs ayia eceree (32) 


where the terms in curly brackets include all those occurring in a except 
SY 
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that containing (¢°,,),,. ‘The difference of these two terms is a single expression, 
that involving G,,,, as in eqn. (23), and differs from that expression only in having 
the negative sign. Hence 


E ace ak cle Bat Gs (33) 


and generally be ‘ = -|. “A seid (34) 
h->0 ee i 


‘These descriptions suffice to determine all the lattice sums we shall need. Since 
the lattice under discussion is symmetrical in the positive and negative charges, 
and as there are only nearest neighbour repulsive forces, it follows that 


E dius E sale sige ons (35) 


where y, 6 take any of the values x, y, z. It can further be shown from the form 


of the lattice sums that 
= : FP aon hee 36 
is aKa E h>0 ( ) 


$6. SPECIAL WAVE DIRECTIONS 


The main purpose of this paper is to evaluate the ‘limiting’ frequencies 
associated with long waves travelling in different directions in the crystal. To 
specify these directions we fix certain relations between the values of h,, h, and hs. 
These directions are then defined in terms of the vectors of the basic cell, but the 
direction of the wave vector s can be expressed without difficulty in terms of the 


orthogonal axes: 

s_.(-h% sa gn yiieemal aE) ; +h, —hg) 37) 

A i; 2a 2b I, Ic = 
where A is the wavelength and the unit vectors i are, as in §2, taken parallel to 
the orthogonal axes of the orthorhombic cell. 

The three cases which will be considered are chosen with a view to simple 
solutions of the determinantal equation for the frequencies (which follow 
from (6)) and are in simple directions in reciprocal space. 

(a) 1. hy=h,=h->0; hg =0. 
2. hy =hz=h>0; h,=0. 
3. hz=h, =h->0; hh, =0. 
These represent waves travelling along the orthorhombic axes in reciprocal 
space, and hence also in the crystal. 
(b) 1. hg= —h,=h>0; hg =0. 
2. hg= —hy=h>0; h,=0. 
3. hy = —h,=h>0; h.=0. 
Here the waves are directed along the face diagonals of the orthorhombic 
reciprocal lattice cell. This does not coincide with the corresponding direction 
in the real cell. 
(c) 1. h,=h,=h,=h-0. Sulie= i, =O en, =h-0, 
2. hy =h, =U; ha =h0: 4 ho=—h,=0; hi =h-=0. 


PROC. PHYS. SOC, LXV, 6—A 27 


402 M. Blackman 


Here the directions are parallel to the (four) space diagonals in the reciprocal 
lattice, which again differ from the corresponding directions in the crystal. 

The solution of the equation for the frequencies is particularly simple in 
case (a), and only in case (c) is it necessary to solve a third-order equation. These 
latter sets of solutions are all identical—which can also be seen from the 
symmetry of the crystal. The values of « and f are chosen to be the same as for 
thallium fluoride, i.e. «= 1-13, 8 =1-38 and f is taken as 1-88. ‘The solutions are 
expressed, in units of B =$7e?/2a°, as jw, where p is the reduced mass and w is 
the angular frequency. The numerical solutions are given in the table. 


(a) Mee) (c) 
Case 1 2 3 1 2 3 
po? 405 3-42 2-96 3574 B10" 5 9954 3-64 
po? 1:03, 1-65" 4-03 0:55 Se 08 At 1-32 
pws” O56) 70:56... ac65 1°32 A205 ear Be 0-68 


The above values fall into two groups, one ranging from 2:96B to 4:05B and 
the second from 0-558 to 1-65B—the first group consisting of the longitudinal 
vibrations. Inserting the value of (the reduced mass of thallium fluoride), 
the frequency ranges are respectively 4-8 x 10! to 8-2 x 10" and 10-9 x 10™ to 
E29 Oe: 

The above set of frequencies forms a sufficient demonstration of the main 
point of this paper—that the ‘limiting’ frequencies are dependent in general on 
the direction of the wave. Even when no distinction is made in the repulsive 
forces between bonds along the different axes, the spread in the frequency of the 
longitudinal vibrations, for instance, is about 15%. 


§7. RELATION TO INFRA-=RED FREQUENCIES 


As has been pointed out in the introduction, the above calculations apply 
insofar as only Coulomb interactions need be considered, and do not necessarily 
include those vibrations which are associated with electromagnetic waves. The 
particular vibrations of the lattice which are ‘ optically active’ (in the sense of being 
associated with electromagnetic waves inside, and outside, the crystal) have been 
considered by Born (1923, p. 728) who has given a general method of calculating 
the electromagnetic lattice sums. In the particular case of crystals of the rocksalt 
type, these optical frequencies coincide with the ‘limiting’ frequency for 
transverse waves (cf. Kellermann 1940 and Born 1923). It might therefore be 
expected that the corresponding calculation for the orthorhombic lattice would 
lead to (three) frequencies which would coincide with certain of the frequencies 
listed in the above table. For the orthorhombic crystal one would further expect 
these limiting frequencies to be associated with those directions which lie along 
the principal axes of the dielectric ellipsoid, i.e. along the orthogonal axes of the 
crystal (case (a), § 6). 

Since the lattice sums used by Born are calculated using the formulation 
of Ewald, these lattice sums can also be used as a check on the accuracy of the 
numerical calculations used in the table. The electromagnetic lattice sums are 


given by 
ae _ €n€x [Op 
xy > | DxO yA aia heal emene (38) 
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where v is the cell volume, €,, the ionic charge and yy =, +H, where 


4 Lf2 
hy= Zw exp | - HEE iq nh, 
y Ceol =r |) 7 
fo = eer os ero Prt (39) 


Here G(x) =1— F(x), where F is the error function. ‘The lattice vector is r’ and 
the reciprocal lattice vector q'. The Ewald separating parameter ¢, is taken as a} 
in the calculation. In the summation for %, the zero term is to be omitted. 
With the same notation and lattice constants as before, the calculation for the 
orthorhombic lattice leads to 


127EM 127M EM 
| ] dank [ | Seas cal 020 ae (40) 
a yy [| zz 


Using the same repulsive forces as before, the three optical (angular) frequencies 
follow immediately (cf. Born 1923) as 


: ak bale 
pwr= an : 
YY Ah=0 Be 


‘These are ww,?=0-550B; pos? =1-05B; ww,?=1-66B, which agree within the 
error of calculation with the transverse frequencies listed in the table, case (a). 


S38. OTHER CRYSTAL TYPES 

A number of other lattices have been considered. Cubic ionic lattices of the 
calcium chloride and zinc blende types yield, in the long wave limit, two transverse 
waves of equal frequency (v,) and one longitudinal wave of frequency 1, as in the 
case of rocksalt. The relation between the frequencies of these waves is the 
same as that deduced by Mott and Frohlich (1939) from electromagnetic 
considerations: ae Ai 2i age (41) 
where v is the volume of the rhombohedral cell, u is the reduced mass, and €« is 
the ionic charge. 

In calcium fluoride we have a more interesting case, as there are three particles 
in the rhombohedral cell and the determinant from which the frequencies are 
to be obtained is of the ninth order. The solutions for wavelengths long 
compared with the size of the crystal have been given by Born (1923, p. 627) 
as two triply degenerate vibrations, one of which is optically active. Calculations 
were carried out using a power law for the repulsive forces, and the results were 
formulated in the terminology used by Born. The following limiting (angular) 
frequencies were obtained: 


2 1 
w,27=v0(D—A’) (— + -) (twice) 


292A’ +D) Gr - ~) (once) aml) \aesser (42) 


) 


w,7=0v(D+2D' le — - (three times) 


where A’ = 8z7e?/3v?, m, and mz, are ad masses of calcium oe fluorine ions, 

D and D’ are the lattice sums defined by Born, v is the volume of the rhombohedral 

cell. The first of these frequencies is associated with transverse waves and is 
27-2 
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the same as that found by Born, the third is also that given by Born, the wave 
being transverse and having zero electric moment per cell. The second 
frequency is associated with a longitudinal wave (with a definite electric moment 
in the reference cell) and is analogous to the longitudinal wave in the rocksalt 
crystal. Terming #,=«, and w,=«,, it follows from (42) that 


2 

eee (= 4 =) pra (43) 
o \m, Ms 

This is a relation of the same form as (41). Applying the method of Mott and 

Frohlich (1939) to calcium fluoride, the relation (43) is also found. 

Finally, the zinc oxide lattice was examined, as an example of a group of 
crystals with ionic forces which do not possess cubic symmetry. ‘The simplest 
cell contains four molecules, and hence the determinant from which the 
frequencies are to be obtained is of twelfth order. In the long wave limit, and 
choosing suitable directions in the crystal, the determinant can be split up so that 
one does not have to deal with more than third-order determinants. The results 
of the calculation confirmed those obtained for the orthorhombic lattice, namely, 
that the limiting frequencies are not unique. In view of this and of the rather 
complicated nature of the calculations, it does not seem desirable to record them 
in detail. 

§9. CONCLUSION 

The general conclusion to be drawn from the above calculation for the 
orthorhombic lattice has already been indicated in the introduction The fact 
that the limiting frequency, in a particular frequency branch, is not unique means 
that the surfaces of constant frequency are strongly distorted in this region. 
The consequence is that where a sharp peak might have been expected in the 
vibrational spectrum by analogy with the rocksalt case, this peak will be 
extensively broadened. As regards the longitudinal branch, characterized by a 
long tail at its high-frequency end, this distortion of the surfaces of constant 
frequency will mean a more extensive and tenuous tail. 

The general features discussed above are expected to apply to all crystals 
which do not possess cubic symmetry. 
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ABSTRACT. Bymaking use of the Poisson—Boltzmann equation an expression is obtained 
for the free energy of the electric double layers of two parallel plates immersed in a large 
volume of an electrolyte of the 1-2 valency type. Two cases arise depending on whether 
or not the divalent ion and the surface charge have the same sign. ‘The formula for the 
free energy is expressed in terms of elliptic integrals. An approximate form for the free 
energy, which is more convenient for numerical computations, is also obtained. This is 
particularly suitable at large surface potentials and is derived for a symmetrical electrolyte as 
well as for the asymmetrical ones. Numerical tables of the free energy as a function of 
plate separation are given and some applications are briefly indicated. 


§1. INTRODUCTION 

DERIVATION of the expression for the free energy of the electrical double 
Ase of two parallel plates immersed in a large volume of an electrolyte 

of the binary, symmetrical type has been given by Verwey and Overbeek 
(1948) in their recent book. Their discussion has been restricted to the case of 
symmetrical electrolytes owing to the mathematical difficulties involved in the 
asymmetrical case. In the present paper we obtain corresponding expressions 
for the free energy of the double layers of two parallel plates in the two cases of a 
1-2 valency type of electrolyte, namely Case I, when the divalent ion has the 
opposite charge to the surface ion, and Case II, when the charge on the divalent ion 
has the same sign as the surface charge. ‘The solution of the Poisson—Boltzmann 
equation for an electrolyte containing both monovalent and divalent ions has been 
already given by Wood and Robinson (1946) and we shall use a similar procedure 
for evaluating the free energy. Robinson (1948, 1949) has shown that if the 
electrolyte is of the asymmetrical type, then we can express our results in terms of 
elliptic integrals only if there are present univalent and divalent ions (in any 
proportion). However, hyperelliptic integrals are encountered in the case of 
higher valency types of electrolytes. 

The exact formulae for the free energy which are derived in this paper are 
rather involved and their numerical computation proves to be laborious. Conse- 
quently it is desirable to seek various types of expansions which are simpler in 
form and, at the same time, reasonably accurate. Onesuch approximate expansion 
is suggested by table XI in the book by Verwey and Overbeek (1948, p. 82), 
which gives the values of the interaction (free) energy of two plates in a symmetrical 
electrolyte. It is observed that this energy very nearly depends on the value of the 
electrical potential at the median plane only. We shall obtain a series for the 
interaction energy, the first term of which is a function of the median potential ; the 
higher terms depend on both the median potential and surface potential and are 
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quite small, unless either the potential or the separation is small. ‘This series is 
given by eqn. (20) for the symmetrical case, by (35) for Case I, and by (46) for 
CaseII. Thus the main part of the energy depends on one parameter only, whereas 
the exact formula for the energy involves two parameters. Consequently, the first 
two or three terms of the above series are much more convenient for numerical 
computation than the exact expression and are sufficiently accurate, except at 
small potentials or separations, particularly in Case I. Verwey and Overbeek do 
refer to the above property of the interaction energy, but there is no mention of the 
expansion which we derive both for the symmetrical and 1-2 valency types of 
electrolytes. In a future paper (to be referred to for convenience as Part IT) we 
shall develop different expansions which correspond to those obtained in previous 
papers (Levine and Suddaby 1951 a, b) for the case of a symmetrical electrolyte. 

In this paper (and also in Part II) we restrict ourselves chiefly to the mathe- 
matical proofs and the computations of numerical tables of the free energy. 
Applications will be dealt with in a still later paper, although some discussion on 
the major role of the coagulating ion on stability properties of colloids is given. 


§2. GENERAL EXPRESSION FOR FREE ENERGY 


Each of the parallel plates is assumed to have unit area, and the edge effects and 
the outer faces of the plates are ignored. We introduce the following notation : 


R_ distance between the plates. 
x distance measured in the normal direction from the median plane. 
y% potential at position x in the overlapping diffuse layers. 
y potential at the surface of each plate (x= +R). 
Py potential at the median plane (x =0). 
n, volume density of ions of species 7 in the bulk of the electrolyte solution 
(=15 2). 
v, valency of anz ion. 
D dielectric constant of the dispersion medium. 
T temperature. k_ Boltzmann’s constant. e electronic charge. 
p= ey/RT, No = eby/RT, Nm = Ct, /RT. 
y =exp(—%), Yo=exp(— 1), k=exp(—mMm). 
«x Debye—Huckel parameter given by «?=47e2B/DkT, B= n,v2, where we 
sum over z= 1, 2. 
o surface charge density. 
E electrostatic energy per plate associated with the double layer. 
Ap difference between the (ideal) osmotic pressure at any position in the 
diffuse layer and the corresponding pressure in the interior of the solution. 
2F (electrical) free energy of the double layers of the two plates. 
2F, free energy of the double layers when plates are at infinite separation. 
V(R)=2F—2F,, interaction or mutual free energy of the plates, i.e. the change 
in free energy when plates are brought to position R from 
infinite separation. 
A(k, ¢) =(1—R*sin*6)1?,  R(R, z) =[(1 — 2?)(1 — R222) p22, 
F(k )=[ ren - [ees x=sin 4, elliptic integral of the first ki 
, ua wae , elliptic integral of the first kind. 


BR, #)= | (kab = | : 


R(k, tial 
os dz, elliptic integral of the second kind. 
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The form of the electrical free energy employed is derived from the Lippmann 
equation, 
2F=-2{" o dip= —2E-2[ Ap ie eis teed eae (1) 

0 ) 


The second relation in (1) follows immediately from a more general formula derived 
by Levine (1951) which applies to a single colloidal particle of arbitrary shape. 
An equivalent equation was derived by Verwey and Overbeek (1948, pp. 78-81) in 
the case of a symmetrical electrolyte. 

The Poisson—Boltzmann equation for + reads 


dys 47e . 
at aa Sa Engg" te re RS eee oe stare (2) 
a first integral which is () - st, Vm Ro Sale | rae (3) 
since djs/dx=0 atx=0. It is convenient to express (3) in the form 
1 dy y dy ot yes 5 
oe gen 49 | pBM0"—F9 |= 4700), (4) 


where p is so choeen that the coefficient of the highest power of yin Q?(y) is unity. 
The positive root is taken if <0, the negative ifo>0. Then the solution of (2) is 


1 Yo dy 
ES 2 al Sue 5 
2 £00)" (5) 
The electrostatic energy per pais is written as 

gid (es Dm db, | De (kT 
ae Gh ae 8x ela = a(S =) pl, ERTS (6) 
where pe | Pk Feo ae TET EAE Re 7) 
kt y ) 


which follows from (4) on transforming from % to y. The corresponding expression 
for the second energy term on the right of (1) is 


R/2 
~2kT&n, { (ielidee fe eee (8) 
0 


Substituting (3) into the first integral expression for E in (6), it immediately 
follows that (1) can be expressed as 


2F=—4E—-RTRin(R%*—-1). hen (9) 
§3. SYMMETRICAL ELECTROLYTE 
In this case we put the valencies of the ion types v, = — v, =v(>0), the volume 
densities n, =n,=n, and chooseo>0. ‘Then 
Bee 2S ae Ort C21 DR It, etiam a Dee casi: (10) 
OO wineries =n Was ie ame NS (11) 
and (5) yields FeR=R1E(R*, Ar) —F(R®,o)}, ks wee (12) 
where £=sina =(Vo/k)77. if ear ae (13) 
The free energy (9) is 
4nkT 


2F=— 


m= | JeRGH* sh) 5, O00) — 4h (E(R®, 377) — E(R’, ah], 
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and when R-> 00 this becomes 
8nkT 


K 


The solution (12) was originally given by Langmuir (1938). The electrostatic 
energy E was first evaluated by Corkill and Rosenhead (1939) and the free energy 
expression (14) was derived by Verwey and Overbeek (1948, pp. 78-81). 

It is possible to express (12) and (14) somewhat differently. In the following 
sections corresponding alternative forms will also be obtained for the 1-2 electrolyte 
cases and the method of arriving at these results will be essentially the same as for 
the symmetrical case. This method is based on the substitution 


1—<2? 


betay ere ek ee (15) 


t= = pete z=sind, z’=sin¢’, AOICAES 4 (16) 

which leads to the relations 
i(k, a7) — F(R, $) = F(R, $'), auclevaie'e (17) 
E(k, $7) — E(k, ¢) = E(k, 6’) — Rk? sind’ cos g’/A(Rk, p’).  ....+. (18) 


These can be derived from the properties of elliptic integrals (Jahnke and Emde 
1945). In the symmetrical case we replace k by k® and ¢ by « Then (12) 


pecomies jcR=htt-R (Relat) © site =<scose/ (Be a) al aaa (12.2) 
and it is readily verified that the second and third terms in the brackets in (14) may 
be replaced by 
1 + (Ryo) } Q (yo) 
24 ee ARP ER’). eae l4a 
hecorian 2 we 


To obtain the convenient expansion referred to in the introduction, the form (14) 
for 2F is expanded in powers of €. Firstly, the Taylor series for F(R”, «) in powers 
of € is readily obtained, and then substitution into (12) leads to the relation 


2KR =k?" F(R®, dcr) —E —9(1 + Rh?) ESB — 263 (1 +R) + fR* JEP 4 LL]... (19) 
The corresponding series for E(k”, «) and O(yo) in powers of € are also derived 
and, making use of (14), (15) and (19), we obtain 


ve folk) + Bk MCL — RYE Hee Od = he eet ee) 


K 
where —_fy(k) = RY? (2 — h° —3k-2°) F(R®, br) + 4&2 E(R®, dt) — 4.0. (21) 


If k is fixed and €+0, which means that 9 ©, the first term in the brackets on the 
right of both (19) and (20) remains, the higher terms becoming negligible. If 
Yo (OF No) is fixed and «R-> co, thenk-1, F(k®, $7)->0o, and again these higher terms 
can be neglected; in this case f,(k)>0, but it does not vanish as rapidly as the 
higher terms. 

An examination of Verwey and Overbeek’s table XI shows that at vy) =4 and 
xR =0-8, for example, the contribution to V(R) from the higher terms is less than 
3%; at vy9=2 and kR=2:3 it is about 2%. We therefore reach the following 
conclusions. ‘The free energy 2F can be written as the sum of three terms, 
namely (i) 2F',, which depends only on yp or 7p, (ii) the part of V(R) which is a 
function of k or 7,, only, and (iii) a correction term which varies with both k and yo 
and which is very small unless y) or «R is small. If, according to Verwey and 
Overbeek, stability properties of colloidal plates depend in the main on the pro- 
perties of 2F' for «R equal to or greater than 2, then the first two terms in the 


VB) = 
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brackets on the right of (20) will usually give sufficient accuracy. ‘T'able 1 gives the 
values of f)(R) and «R=(«R),, when € =0 (7) = 00) for various k®. We have used 
the tables of the complete elliptic integrals of Fletcher (1940), who chooses intervals 
of the modulus (here k”) instead of the usual angle sin-'k®. We have not 
tabulated the higher terms in (19) or (20), but these are readily computed. It is 
noted that the series (19) and (20) converge for all «R>0, 7) >0. 


Table 1 
ke (KR) « Folk) fae (KR) Folk) 

0-0025 0-3142 27°573 0-4 4-1489 0:23455 
0-005 00-4443 18-436 0:45 4-4566 0:17450 
0-01 0-6283 12-021 0:5 4-7680 0:12920 
0-025 0:9936 6-427 0-55 5-0886 0:09468 
0:05 1-4058 3°715 0-6 5:4245 0:06824 
0-075 1-7231 2°573 0°65 5-7837 0-04801 
0-1 1-9919 1-9259 0-7 6:1769 0:03261 
0-15 2:4473 1-2103 0:75 6:6199 0:02106 
0-2 2°8387 0-8243 0°8 7:1386 0:01259 
0:25 3-1925 0:5860 0-85 7°7813 0:006646 
0-3 3°5231 0-4270 0-9 8-6545 0:002784 
0-35 3°8396 0:3154 0-95 10-099 0:0006583 


§4. CASE I. DIVALENT ION OF OPPOSITE CHARGE TO SURFACE 

We assume the surface charge density o is positive and consider a 1-2 
electrolyte with a negative divalent ion, for example K,SO,. Taking K* as 
species 1, v, = 1 and SO,’ as species 2, v, =2, n,=2n,. Then 


B= On) p= 2400, 67/DRTs p= (2/3) ate (22) 
Oly) > (x — hy" + 3) = {(y—a)(y— Db) ye) ena (23) 
where 
1h-2 b I 1 1 2 a i 
h=k+ 3k, a=k, ee 8 apy, |=: ab+bce+ca=0, abc= —4..(24) 


Since ¢ is positive, %,,<, and hence y<k<1(yy=k at R=0). The disposition 
of the roots and the range of y are b>1>a=k>y2>y)>0>c. If we introduce 
the substitution 


u=sin D={(y—Ci/(a-—c)pte wm wae (25) 
and put 
Uy =sin Do ={(Yo—e)/(a—e)P”, ky={(a—e)[(b—c)p, se (26) 
then (5) becomes 
$xR = {6/(b—c)}¥? {F(R,, 47) — F(Ry, Do)}, nae wee (27) 


which is the solution of the Poisson—Boltzmann equation. ‘The formula (9) for 
the free energy becomes 
2F=—4E—-RTRn,(2h-3),  —— .sseee (9a) 


and in Appendix I this is evaluated, giving 


| 
Wo SunE | BR(6b—2h—3)+ aly 
0 


~ 34/6 (b —c)¥*{E(k,, $7) ~B(h,)} |, arte (28) 
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where Q(yp) is given by (23). The value of 2/,, the free energy of the two plates 
at infinite separation, is readily found. When Ro, then k-1, E(h,,37)>1, 
E(k,,®,)—sin ®,, the first term on the right of (28) vanishes and the last two 
terms yield 


Piuaage mkt | v3 (2+ =) (1 +2y)1—9 |. ae (29) 
0, 


The same expression may also be obtained directly by putting k = 1 in the expression 
for Q(y) in the integral (7) and evaluating (9). If the surface charge density 
o<0 and the 1-2 electrolyte present has a positive divalent ion, it is clear on 
physical grounds that the free energy (28) remains unchanged. It is readily 
verified that the preceding analysis is reproduced exactly if we replace y and k by 
their monigesl (which is equivalent to changing the signs of o and of the divalent 
ion). 

In the preceding section we made reference to an alternative form for the free 
energy based on the substitution (16) and the use of the relations (17) and (18). 
For the case under consideration this substitution will be written 


UN Geer OY (igs Seat ga ela Bee Nc 
i=in = (a) = R, (5) Gg MER Geret fe (30) 
Then (27) can be expressed as 
6 \12 ; 1 /yy—a\12 
ie zal ret cn / eet oy ea 0 sar 
ic (=) FsOe), san, z (25) oe (27 a) 
and the second and third terms in the brackets in (28) become 
6(642 : 
eee O(9)— 3/6 (b—c)¥? E(Ry, Do’). a. s- (28 a) 
Yo(b—Yo) 


It is also possible, although more laborious, to obtain the above results directly 
from (5) and (9) by employing the substitution (30). 

Finally we derive an expansion for the free energy difference V(R) which 
corresponds to ‘ee Firstly, 


ae =|. RE Riko) = F(R, Q)+ aS [ ss 


0 OY) 
—c\12 » 
nesses (vot Haye? — dott BhtyeF + «..), 
Sea tee (31) 
where sinQ=—c/(a—c). This expansion is obtained if we transform from u 


to y, making use of (25) and (26), and then we expand the integrand in (31) in 
powers of y. Similarly 


1 
E(k, ©) = E(Ry, Q) + (b= oie {by —3V9? + gbhy 8 


—2(b + h)yot + H(BBh2 +1) yb + 0.0}, ceceee (32) 


noting that (6—¥y)/Q(y) replaces the integrand 1/Q(y) in (31). It is now con- 
venient to introduce the substitution (16) where z =sin Q, ice. 


; cos Q a(b—c))\ 12 2a+b\12 
ae IN SEO) ace) i (3) tito pi ty 


Double Layers in a 1-2 Electrolyte 4II 


which follows from (24). Making use of (17) and (31), (27) becomes 


F(R yan) y hy 3 y 4 3 
5) R=/ : we CA ec au 0 80" 
Le V ‘6 Ress /2 6x72 =a a iigwoes rk eeeee nO (34) 


To obtain V(R), the quantities O( yo) and 2F.,, defined by (23) and (29) respectively, 
are expanded in powers of y) and use is made of the substitution (33) and of the 
relations (17), (18) and (32). Then it follows from ed (29) and (34) that 


: — 1 
V(R) = E (k) + Z73 th 3) ay 84+ Vo es $)?(4h +3) + re 
_ nt 


B23), ay) = ae ee She Sarre ae (35) 


where (Rk) = (ge) (2h+3—656) F(R,, Q’) 


4 34/6(b cy BCR OV 3 1/6 l Set) Oe (36) 


The equivalent form for f, (A) in terms of Q is readily obtained. 

Just as in the symmetrical case, 7,(R) is the main contribution to g,(k, Vo) 
provided «R or 7p is not too small. Indeed, the expansions (31), (34) and (35) 
converge only if yy<|c|, and this condition implies that «R must be greater than 
0-71, 0-17, 0-07, 0-03, 0-01 and 0-004 approximately when 7 =1, 2, 3, 4, 5 and 6 
respectively. ‘Tables 2 and 3 give the values of «R and the function g,(k, yo) for 
the above range of 7). When y)=2 and«R>1, the first three terms in the expan- 
sion (35) yield V(R) with an error of less than $°%. As either «R or no increases, 
the accuracy of these terms improves; thus we have at least four-figure accuracy 
when 7) =2, «R>3, or np =3, kR>1 or ny =4,neR>}. When no=1,«R=1 the 
error is about 5°% and, although this diminishes as «R increases, the third term in 
(35) is larger than the second, which means that we require the higher terms in the 
expansion. ‘The corresponding first five terms in the series (34) for «R are even 
more accurate; thus when 7) =1,«R=1 we have an error of less than 3°%, and when 

=2,«R>1 at least four-figure accuracy. It is convenient to divide table 2 into 
two sections; in the upper section we need not tabulate g,(k, v9) since it is identical 
with g,(k, 00) =f, (k) if we retain only four figures. In table 3 the values of g,(k, Vo) 
are calculated for small «R. We used the exact formulae (27) and (28) for both 
No = 1 and small «R and also as a check in the range where the expansions (34) and 
(35) are applicable. Where the values of & are not given in table 3, we have applied 
another series for g,(k, Vy), which is suitable at small separations and which will be 
derived in Part II. Also in the case yn) =1 and large «R still another expansion for 
£,(k, yy) has been used; this will also be given in Part II. 

It is observed that the correction to the limiting function f,(k) is of order 
yo? = exp (— 39), which is very small for 7) equal to or greater than 2. This is to 
be compared with the corresponding correction in the case of a 2—2 valency type of 
electrolyte, given in (20), where we find the same factor in €° = k-* exp ( — 3) since 
v=2. The range of «R and 7, over which the series (35) can be applied is much 
greater than that of the expansions developed in Part I]. However, these other 
series can be generalized to higher valency electrolyte types, whereas the corre- 
sponding extension of (35) would involve hyperelliptic integrals. Also, to obtain 
F(R) for large «R, we subtract two quantities which are very nearly equal and this 
requires a double interpolation of the tabulated values of elliptic integrals, which is 
found to be rather laborious. 
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Table 3 
No 
k 0:00248 —_ — 0:00048 0-0095 0-0185 0-036 
6 £i(R, Wo) 694-96 615-0 542°, 282°1 139-6 69-81 34-13 
KR 0 0-002 0-004 0:01710 0:04301 0:09206 0-:1873 
k 0:00674 _ — 0-011 0:0185 0-025 0-03 
5 2(R, Yo) 253-27 226-2 201-2 121-9 70:01 50-64 41-69 
KR 0 0-005 0-01 0:03473 0:07678 0:1124 0-1399 
k 0:0183 — oe 0:032 0-050 0:057 0-075 
4 i(k, yo) 90-810 79-84 69:9, 39:36 23-69 20-33 14-64 
KR 0) 0-015 0-03 0-1066 0:2071 0-2456 0:3441 
k 0:0498 — — 0:07 0-1 0-15 
3 SARS V5) So L-113 27:18 23:60 16°44 10-28 5-850 
KR 0 0-04 0-08 0-1891 0-3638 0-6412 
k 0-1353 — -- 0:17 
» 2i(R, Vo) 9°3315 8-055 6:883 5-148 
KR 0) 0-1 0:2 0-3836 
k 0-3679 — — 0:75 — — 
1 £i(R, Vo) 1:°7669 1-5133 1:2747 0:06968 0-03638 0-01341 
KR 0) 0-2 0-4 3:348 + 5 


§5. CASE II. DIVALENT ION OF SAME CHARGE AS SURFACE 
It is convenient to.consider again a 1—2 electrolyte with a negative divalent ion, 
but this time to take the surface charge as negative. Since the potential ¢ is 
negative, the distribution of the roots a, 6, c and the range of y are now 
Vo zy za=k>1>b>0>c. From physical considerations one would expect 
Case II to follow directly from Case I by replacing 7) by —7. ‘The details of this 
method of treating Case II are as follows. We introduce the (real) substitution 


: a—c\i2 J 
vU=sin p= (=) = are arelsiiallene (37) 
and put 
a—c\v2 b—c\12 1 
Vp =sin By = (<=) ke = (=) = 7 Parnes. (38) 


Then, making use of the first equation in (31), we can write 


F(k,, 4) — F(R,, ®o) =e iF = ko F(R, $7) — F(Re, Bo)]. -. ++. (39) 


Thus (27) becomes 
£cR-= {6/(a—c)}¥? {F(Ro, 47) — F(Ro, Bo}. = we ee - (40) 
This can also be obtained by introducing the transformation (37) directly into (5). 


In equation (27) ®, is real and k, <1 if the potential is positive, whereas in equation 
(40) By is real and k, <1 if - potential is negative. Similarly 


El. ae v) 1 (7? A (Ro, B) 


ie eal 
a, sin?e- 


o(1 — v2) hy. 


dv 
R(Rp, v) 


E(k, 377) - Ey, ®,) ae dp 


= 7M Bo) cot By + 7, (El 377) — E(Re, Bo)} 


+ OS lke, bn) Fly Bh ee (41) 
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where the final form is obtained after integrating by parts. Using eqns. (24)—(27), 
it can be verified that (28) becomes 


a «R(6k —3 — 2h) + (2+ £) O10) 


= 3/6 (a~€)!* {Bl $n) — Elbe Bo)} | ---(42) 


where Q(y) is again given by (23). This is the required formula for the free 
energy in Case II. We find 2F,, in the same way as in Case | and obtain a form 
identical with (29). It is noted, however, that in Case I y)>7,>0 and thus 
Yo<k<1, whereas in Case II 49<y,<0 and thus yy»>k>1.‘ In Appendix II 
it is shown that we can arrive at the result (42) by starting with the original formulae 
(9) and employing a method analogous to that used in Case I. 

The alternative expression for the free energy is readily obtained. We 
substitute 


Pa 


yea, J —y? \1/2 » — q\ U2 
O-=— Sin B = (eo) = (2 = ;) + By svexsisheke (43) 
Then (40) reads ; 
6 \12 Yo —@\12 
ak oo (eae / a 0 
LR (; = .) F (ky, Bo), sin Bo & =) ee 2 (44) 
and the second and third terms in the brackets in (42) are transformed to 
/O(2yy +b 
Vv y 
VOC¥0* 9) (34) —3{6(a—c)}"* Blas Bo") eae (422) 
Yo(¥o — 5) 


which corresponds to (28a). As in Case I, this equivalent form for the free 
energy has also been derived by the more laborious method of introducing the 
substitution (43) directly into (5) and (9). 

Finally, using (37) and (38), we expand F(R, By) in powers of 1/y» and then (40) 
becomes 
jy 8 Lie ha) | 1| * a 

* J O() 
PY, ky, 377) Dh 3 ee 1 (5h?—4) 
6 =v6| 2 (a—c)¥2 = 6y 32 40 y92? wie Yor? came Ft) 

which tae to (34). Similarly, we may expand E(k,, By) in powers of 1/79, 
noting that in place of 1/Q(y) the integrand is now (y—)b)/(y—c) O(y). Making 
use of the latter expansion and of the corresponding ones for O(y)o/(vg—c) and 
ZF, we derive from (29), (42) and (45) that 


4n.RTT . 1 1 1 
V(R)= — EG ae eee cE alle 
(R) fall) + seh coke pth 9) (4+ 3) as 
we 1 4n,RT 
+ 3 A — 3)? (4h? + 4h 3) yi +. = —— 82(k, Yo): say, (46) 
where 


fx(k) = Cres 5) (2h + 28-43) Fk, bn) + 34/6 (ac) Bll, x) —9. .. (47) 


Here again f,(k) is the main contribution to g,(k,y) provided «R or no 1S 
not too small. The series (45) and (46) converge for all «R>0, 1 >9 since 
|c|<b<a<yo. However, the correction to the limiting function f,(k) is now of 
order exp (—3 7/2), so that for large 7, the series (46) does not converge as rapidly 
as the corresponding one in Case I. We note that the expansion (46) is to be 
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Table 5 
No k 403-4 —- -—- 200 100 90 
6 B(R, Vo) 89:582 79°67 70:49 47-26 30°34 26:23 
KR 0 0:05 0-1 0:2735 0:-5142 0:°6075 
k 148-4 _- —— 100 80 
5 (Rk, yo) 50-984 43-78 BO 35 32-01 DAS 
KR 0) 0-1 0-2 0:2978 0:4088 
k 54-60 — — —- 
4 gk, yo) 27-697 25-06 22°54 20-22 
KR 0 0-1 0-2 0-2 
k 20-09 — — 
ao edhe is-78l, 11-94 10:23 
KR 0) 0:2 0-4 
k 7-389 © — — = 
» (Rk, Vo) 5:°6911 5-106 4-538 4-002 
KR 0 0:2 0-4 0-6 
k 2:718 — — — 2°65 2°6 2°25 
1 2(R, Vo) 1:4051 1-3409 1-2769 41-2135 1-0940 0-9952 0-5869 
KR 0) 0-1 0-2 0-3 0:-4922 0:6570 1-458 


compared with the series (20) for the case of a 1-1 valency type of electrolyte. In 
tables 4 and 5 we have calculated «R and g2(k, yo) for various 77. When the values 
of go(k, yo) are not shown it is understood that g.(k, yo) =f2(k). The rapidity of 
convergence of the series (45) and (46) follows a similar pattern to Case I, the differ- 
ence: being, of course, that to obtain the same accuracy for a specified y) and the 
same number of terms we need to choose a larger«R. For example, if we retain 
the first three terms in (46), the error in V(R) is less than 4% when 7 =2, 
KR >2'5 and about)», if 7, = andxkR=2Z. 
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Verwey and Overbeek (1948, p. 33) demonstrated the small effect of the ion 
carrying the same charge as the surface on the capacity of a single flat double layer. 
In the case of two parallel plates, however, their calculation only shows that-this 
ion has a minor effect on the free energy at infinite separation, or on the repulsive 
potential V(0)= —2F,, at contact. The influence on V(R) at finite R, and parti- 
cularly at large R, remains to be evaluated. As the above authors point out, the 
colloid chemical facts show the secondary importance of the ion of the same charge 
as the particle. ‘This property can be investigated theoretically by applying the 
results of the present paper, and some tentative conclusions will be described here. 

Let the interaction energy for 2-1 and 2-2 valency types of electrolytes be 
denoted by V2,(R) and V..(R), respectively, in the case of a divalent coagulating ion 
(ie. the ion carrying the opposite charge to the plates); the corresponding 
Debye—Huckel parameters are «(2,1) and «(2,2). Assuming the same plate 
potential ys) and the same concentration m, of the divalent coagulating ion (so that 
(2, 2) =(4/3)«?(Z, 1)) in these two electrolyte types, we have plotted V,,(R) and 
V..(R) as functions of R. It is found that, except for large R, the two curves 
obtained do lie quite close together, particularly for large y)._ Thus at R=0 the 
ratio V;,(0)/V2,(0) has the approximate values 1-06, 1-024 and 1-01 at ms Dig a: 
respectively, and tends to unity as 7 increases. As the plates are separated the 
ratio Vy(R)/V.,(R) initially increases somewhat, but then the above two curves 
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cross each other, and for large R («(2,2)R<3, say) we have the approximate 
een Vo(R)/Vox(R) = A(mp) exp [—0:134(2,2)R], eevee (48) 
where A() increases from 1-48 to 1-79 as 9 varies from 1 to 6. ‘The relation (48) 
yields an upper bound for the value of «(2, 2)R at the crossing point where 
Voo(R) = Voi (R). 

In the case of a univalent coagulating ion we denote the energy functions for 
1-2 and 1-1 electrolytes by V;,(R) and V;,(R) and the parameters by «(1, 2) and 
x(1,1) respectively. Assuming again the same plate potential 4), but now the 
same concentration 7, of the univalent coagulating ion (so that x? (1, 2) =1-5«(1, 1)), 
we have plotted V,.(R) and V,,(R) as functions of R. The two curves are more 
separated than in the previous case, although still never differing very much 
except at large R. At R=0, V,,(0)/V,,(0) has the values 1-126, 1-07, 1-04, 
1-025, 1-015 and 1-01 at y)=1, 2, 3, 4, 5 and 6 respectively, and the behaviour of 
V,.(R)/V,,(R) as a function of R is similar to that of V.9(R)/V.,(R). Thus at large 
separations («(1, 1)R>4, say) we have approximately 


Vi2(R)/Vii(R) = Blo) exp[—0-225e(1,1)R], «see (49) 
where B(n 9) increases from 1:63 to 2-64 in the range 1....6 for m9. Again (49) 
provides an upper bound for the value of « (1, 1)R at the crossing point. 
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In obtaining (28) we simplify the notation by considering the indefinite 
integral, J say, which corresponds to the definite integral J, as defined by (7). 


Then J can be written as path py PTGS I eaves e 
: ydy Pig aaa oe, rt 
nee far eo) ~ (b—c)P | A, 2(b—c) | Gd Dees (A 2) 
“d® 
=l5¢ O(y) = aoe) & Pe Ae Poe (A3) 


dy 2 ; du 
a is FOG)" ae=-o"| (rm ae 
A, =A(R1,®), m=(a—c)/c. 
J.,is not yet in a standard form but, adopting the usual procedure, we differentiate 
the expression wR(k,, u)/(1+ mu?) and then integrate to obtain a reduction formula 
which expresses J, in terms of the three kinds of elliptic integrals. Writing the 
result in terms of the variable ®, the integral in (A 4) reads 


; d® 1 (m+k,?)pd® 1/7 
ee ee aS ge fii) 
| cement OR, mal 2 | Ieee A ee 
' d® A, sin ® cos ® 
+¢| (1+msin?®) A, % Aer | ieee c3) 
; ki? 1+hk,? 3k, 
where ee +, C=1+ sal ee (A 6) 
m m 2m m 
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If we substitute for k, and m we find that ab+bc+ca=0 appears as a factor in C 
and hence C=0; also./, simplifies to ; 
J 2b d® 

EE SLE hoi 
The integral (7) for J, and hence the expression (6) for £, are immediately obtained 
by introducing the limits of integration $7 >® >, when k >y > yy in the formulae 
for J,,J,andJ3. It is convenient to eliminate the elliptic integral of the first kind 
by means of (27), and then it is readily verified that the formula (9 a) for the free 
energy becomes the required expression (28). 


2(b — 2 i A,d® +2 a eee (A7) 
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To determine the integral (7) for J in Case II we _ substitute 
#2 =1/(y—c) =v?/(a—c) and write 
S(t) =[1 —(b-—c)#}{1 -—(a—c)t} Pr. (A8) 


Considering the indefinite integral J, we have 


1 
J=~2H,+ 5 ihy+ ais, ee eeee (A9) 
dt dp 
= | SS Saige | — 
where H, | ESO) (a—c) le Sages See ey glace (A 10) 


= —(a—c)!? E cot B— (= a [a.ae | » Ay =A (hy, B), 


dt s 1 dp 
i= | (l+cf?)S(t)  (a- aa (ht singB) AS? Aten seks 
which is an elliptic integral of the third kind, and 
z ds aia dp : 
Ham | reas ~ aot | (eaten FS «(AT 


This can be determined by applying (A 5), and we now find that 
D=1+(1+hq?)m + ky? = —c?, C=14+2(1 +h?) + 3k,2m? = —2/c3. 
Making use of these results, (A 9) becomes 
b\ dB ; c\ O(y) 
Aes - (2+ 2) wt alg ole Aa hh) ae 
(a—c) ( o) ise (a—c) | dB + (2+ 5) 20), Space (A 13) 
noting that the coefficient of the elliptic integral of the third kind vanishes. Intro- 
ducing the limits of integration $7>8>@®) when k<y<yy into (A 13), we 
readily derive that the formula (9) for 2F leads to the required result (42). 
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ABSTRACT. The rotational analysis of the near ultra-violet band-system, p!II —x!X*, has 
been extended. Twelve bands have been analysed, namely the 7,1, 6,0, 6,1, 5,0, 5,1, 
4,0, 4,1, 3,0, 2,0, 1,0, 1,1 and 0,0 bands. The constants derived for p!II are 

B,, =0-2664,—0-0021, (v’+-4); De [2:91 —0-03(v’+ 4)] x 107", 

G,, =513-1,(@' +4) —2:93(v’ +4)?» ve=35026°8, cm-"; re =2:058, A. 
A number of perturbations have been found in the upper state. They appear to arise from 
interactions with two electronic states, a 9X state and a singlet state, II or 1A. It has 
not been possible to determine the absolute numbering of the vibrational levels of the 
perturbing states, but in all other respects fairly complete information has been obtained 
about the °X state. The perturbations are similar to those already known in the spectrum 
of CO, and the analogy between the electronic states of the molecules is thus 
strengthened. 
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situated in the ultra-violet region of the spectrum; the E-x system, which 

lies in the region 2050-2600A (Vago and Barrow 1946), and the 
D—x system, which extends from about 2500 to 3900A (Barrow and Jevons 
1938). These are both obtained in absorption, and the D—x system has also 
been obtained in emission. Vibrational analyses are given in the papers 
mentioned above. A further system has been obtained in emission in the visible 
region (Barrow and Jevons 1938), but it has not been analysed. A rotational 
analysis of the 0,1, 0,2, 0,3, 0,4 and 1,5 bands of the D—x system (Barrow 
1946) showed that the transition was p1II-x!*. 

The vibrational analysis of the D-x system suggested that the heads of the 
bands with v’=7 are abnormally situated, indicating a perturbation in this 
level. An extension of the rotational analysis of this system to higher vibrational 
levels in p1II was therefore undertaken with the object of getting information 
about this perturbation and at the same time of determining better values of the 
constants. This paper presents the analysis for all the upper levels from 
v’ =0 to 7: it shows that two electronic states perturb the 'II state, one of which 
is a 3” state, the other a singlet state. The °X state perturbs the "II level in 
such a way that the interactions from the substates °X,_, and *, , , are visible 
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TT sist band-systems ascribed to the SiS molecule are known. ‘Iwo are 
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in the R and P branches and the interactions from the substate °X; in the 
Q branches. The very strong perturbation causing the head-displacement of 
the level v’ =7 arises from interactions between the 'II and the *X° states. 


§2. EXPERIMENTAL 


The spectrograms were taken in absorption in a first order of a 21-ft Wood 
concave grating (165000 lines) having a dispersion of about 1:25A/mm. ‘The 
absorption tube was made of quartz, with quartz windows. ‘The length of the 
tube was 75cm, of which 50cm were heated in an electric furnace. ‘The 
temperature of operation was about 1000°c. The tube contained initially 
a mixture of Al,S; powder and quartz sand, which gives SiS on heating. 
A hydrogen lamp was used as source of continuum. The exposure time on 
Ilford Zenith plates was about 30 minutes. ‘The band-lines were compared 
against Fe arc lines whose wavelengths were taken from the M.I.T. tables (1939). 


$3. STRUCTURE OF THE BANDS 


The appearance of the structure of some of the bands is shown in the 
spectrograms reproduced in fig. 1. A Deslandres scheme of band-origins is 
given in table 1. ‘The wave-numbers of the band-lines are given in tables 2—5, 
where overlapped lines are indicated with an asterisk. Perturbations have 
been found in all the upper state vibrational levels examined from v’ =0 to v’ =7. 
Those in the level v’=0 are, however, very weak, and they were not detected 


Table 1. Deslandres Scheme of Bamd-Origins 


7 37596-41 
476°67 
6 | 37864:08 744:23 37119-74 
477°68 47810 
5 | 37386:40 744-76 36641-64 
48386 483°50 
4 36902:54 744°40 36158-14 
489°74 
3 | 36412-80 
49598 
2 | 35916-82 
501°23 
1 | 35415-59 744-45 34671-14 31744-61 
506°95 507°26 


0 3490864 744°76 34163°88 739:34 33424:54 734:23 32690-31 729-04 31961-27 


v'-> 0 1 2 3 4 5 


in the original analysis (Barrow 1946). A strong perturbation in v’ =7 
displaces the origins of bands with v'=7 by about 5cm™. In tables 2-5 
romans and italics are used to distinguish interacting levels—except that the 
type is not altered at the weak perturbations in v’'=0. For the weak 4,1 band 
only the Q branch is given. 

Silicon has three stable isotopes, sulphur four, but of these 28Si (92-2 %) and 
325 (95-1%) are much the most abundant. We have made no effort to measure 
lines other than those attributed to these isotopes. 


J QJ) 
4 34908-81* 
5  08-27* 
6 07-74% 
7  07-01* 
8 06-23* 
9  —05-49* 
10  04-68* 
11 «03-43% 
12 02-54% 
13. 01-47* 
14 —00-24* 
15 $99-08* 
16 —-:-98-03 
17:96 -76 
oS. 95-45 
19 =: 93-96 
20 = 92-48 
21 90-89 
22 ~=—- 89-27 
23 4230 s«887-66* 
24 = 85-80* 
25 83-92% 
26 81-89% 
27. 2S s«79-84* 
28 77-70* 
29 75-56* 
30 73-18* 
31 70-79% 
32-68-51 
33 66-06 
34 . 63-45 
35 60-91 
36 58-15 
SWf 55-39 
ae «52-45 
39 = 49-48 
40 46-56 
41 43-43 
42 40-17 
aa «637-03 
44 33-61 
45 30-28 
46 26-73 
a7 «93-19 
48 19-57 
49 15-88 
50 12-00 
50 
51 08-15 
51 
Be «04-22 
53. 00-18 
54 796-13 
55-91-92 
56 87-71 
Mee 83-32 
58 78-97 
59-74-39 
50 69-82 
61 65-20 
62 «60-54 


63 55-69 
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Table 2. Wave-numbers for the 0,0, 1,0 and 1,1 Bands 


0,0 Band 
R(J) 


34908 -81* 
08-27* 
OT 
07-01* 
06-23* 
05-49* 
04-68* 


03 -43* 
02-54* 
01-47% 
00-24* 
899-08* 
97-50 
95-98 
94-62 
93-03 
91-36% 


SOS 
87-66* 
85-80* 
83-92% 
81-89* 
79-84% 
ATO 
15-56" 
MAB 
FO-79 


68-11* 
65°55* 
62-92% 
60:24* 
57-44* 
54-55* 
51:65" 
48-52* 
45-46* 
42-28% 


39-01* 
35-66% 


B2-257 
28 -83* 
2pe22," 
21-41 
17-66 
13-84 
10-01 
05-94 


02-08 
797-76 
93-60 
89-20 


PY) 


34891 -36* 
89-59% 
87-66* 
85-80* 
83-92* 


81-89* 
79-84* 
HEIs 
TERS 
WA8* 
870-79* 
68-11* 
65°55* 
2-92 
60:24* 


57:44* 
54-55* 
51-65% 
48 -52* 
45-46* 
42-28% 
39-01* 
35-66* 
S22 50 
28 -83* 


PAIL» 
Zits 
18:07 
14-24 
10°38 
06°51 
02-50 
798-33 
94-23 
90-08 


85°73 
81-35 


76°85 
72:16 
67-70 
63-07 
58-16 
53-24 
48-32 
43-22 


58°35 
Soll 
21-93 
22:54 


QV) 


35413-30 
UPAR Tigo 
ilalee}ens 


11-02% 
10-28 

09-39 

08-33 

07-374 
06:09* 
04-74* 
03-42* 
01-96* 
00-36* 


398-86 
O7-12 
95-47 
93-62 
91-70 
89-67 
87:65 
85-45 
83-28 
80-96 


78:62 
76°15 
73-60 
70-95 
68-29 
65°65 
62:88* 
59-76* 
56:69* 
53°54* 


50-30* 
47-03* 
43 -62* 
40-19* 
30-75" 
33-02" 
29-33* 
25-56% 
21-67* 
17:68* 


13-47* 
17-68% 
07-18 
10:56 
05:96 
01-56 
297-14 
92-76 
88-27 
83:66 
78:97 
74°20 


69-38 
64-43% 
59-38 
54:35 


1,0 Band 
RW) 


35415-65 
15-03 
14-51 


PV) 


13°77 s03 9003) 


Ne fife 
(edie 
102% 


09-86 
08-76 
O7-378 
06-09* 
04-74* 
03-42* 
Ol: 96% 
00-36* 
398-54 
96-70 


94-88 
2-97 
90-89 
88-79 
86-60 
84°35 
82-04 
USD 7 
77-06 
74-45 


TACs) 
68-99 
66-11 
62-88* 
60-47 
57-29 
54:07 
50°86 
47-53 
44:07 


40-53 


36-75* 
33-02* 
30-15 
25-97 
22-02 
18-12 
13-98 
09°85 
05*55% 


01-14* 
296-74* 
92:27* 
Sie 


94-48 
92-39 
90-42 


88-35 
86-16 
83-95 
81-52 
79-10 
76-61 
74-04 
71:35 
68°57 
65-65 


62:88* 
59102 
56-69* 
53-54* 
50-30* 
47-03* 
43 -62* 
40-19% 
36°75* 
33-02* 


29-33* 
25:56* 
PH yes 
17-68* 
13-47* 
09-35 
05+55* 
01-14* 
296-74* 
92:27* 


87-72* 


83-06* 
78-34% 
73-49% 
68-49% 
64-43% 
59-09 
53°88 
48-72 
43-41 


38-02 
32-62 
27-09 
z1e53 


QV) 


34665 -09* 
64-08* 


1,1 Band 
RV) 


34668 -14 
67:02 


66°13 
65-09% 
64-08* 
62°72 
61-53 
59°89 
58-56 
56:97 
55075 
53-45* 


Sylezisye 
49.94% 
48-10* 
46-12* 
44-00 
41-68 
39°52 
37-11 
34°77 
32-26 


29-76 
27-09 
24-41 
PWNS) 
18-95 
15:95 
12-83 
09-80 
06-50 
03°25 


699-90 


96-40 
92-63 
90-02 
85:90 
82°17 
78°36 
74°53 
70°46 
66°33 


62:24 
57°92 
53-68 
49-22 


421 


P(J) 


3465407 
52-23 


WoNHAU & 


10 
11 
iW72 
13 
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Table 2 (cont.) 


1,0 Band 
RV) 
83-06* 
78-34% 
73-49% 
68-49% 
63°59 
58:56 


53-3 
48-12 
42-80 
37-45 
31-88 
26:25 
20:64 
14-87 
09-09 
03-06 


197-01 
90-87 
84-81 
78-46 
72:10 
65-72 
59-05 
52-49 
ms )ooyl 
SOro7) 


31-94 
24-97 
17-86 
10-84 
03-50 
096-10 
88-88 
81-28 
73-68 


PJ) 
15.74 
10-05 
04-14 
198-22 
92°21 
86°14 


79-94 
73-71 
67-43 
61-01 
54-40 


47-64 


41-17 
34-32 
27-49 
20-50 


13-49 
06-39 
099-32 
91-97 
84-54 


QV) 
10-93 
05-82 
OORT 9% 
495-57* 
90-26% 
84-85* 


79-36% 
73-66* 
68-13 
62-39 
56-66 
50-78 
44-95 
38-79 
32-67 
26-48 


20-24 
13-97 
07-38 
00-85 
394-24 
87-51 
80-71 
74-09 
66-88 
59-72 


52:96 
45-47 


1,1 Band 
R(J) 
44-79 
40:21 
35°52 
30°88 
26:09 
21-30 


16°37 
TA35 
06-03 
00-79% 
495:57* 
90-26* 
84-85* 
19-308 
73:66* 


Table 3. Wave-numbers for the 2,0, 3,0 and 4,0 Bands 


0,0 Band 
QV) RV) PV) QU) 
50-79 84:91 1712 49-22 
45-74 80:52 11-62 43°93 
40-63 76:02 06-12 38-61 
35°67 71:28 00°31 33-16 
30:34 66:59 694-67 27-63 
25°19 61:97 88-92 22-06 
19-77 57:04 83-03 16-41 
14-24 52-03 77°28 10-60 
08-83 47-16 70:95 04-80 
03-20 41-69 65-09* 198-87 
697-40 36°81 58-56 92-86 
91°72 31-48 52-23 86-76 
85:90 26:16 80-53 
79°85 20°59 74-39 
73-84 15-14 68-03 
67:70 09:61 61-53 
61°53* 04-03 55°05 
55:07* 698-23 48-46 
49-01 92-49 41-70 
42-32* 86-57 34:90 
36-03 28-07 
29°38 21-14 
14-03 
06-89 
099-79 
92-42 
85:12 
77°30 
69-77 
62-15 
54°37 
46°63 
38-60 
30-65 
22-49 
14:15 
05°89 
34997 -56 
89°20 
80°59 
71:98 
63:24 
54:40 
45-45 
36°70 
27-60 
18:04 
2,0 Band 
OO) SREP) SOU) 
35912-07 
11-21 35917-61 
10-24* 16:97 35903-92* 
09:01* 16°37 02°31* 


3,0 Band 
RY) 


P(J) 


QV) 
36901 -25 
00:54 
899-81 
99°25 
98-55* 


97:54* 
96°50* 
OSE 27e 
ayeokes 


4,0 Band 
RV) 


P(J 
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Table 3 (cont.) 


2,0 Band 3,0 Band 4,0 Band 

J Qi) R(J) P(J) Qi) R(J) P(J) QJ) R(J) P(J) . 
14 08-04 15-65 00-83* 92-90 
15 06-70 14-98 898 -94* 91-45 
16 05-33 13-95 97-16* 89-95 36898-55* 
17 03-92* 95:06 36399-41 88-39 97-54% 
18 02-31* 13-02 92-91 97-79 36388 -64* 86-72 96:50* 36877-72* 
19 00-83* 11-62 96-13 86-43* 84-99 95:27* 75-41* 
20 898-94* 10-24% 89-54 94-25 84-20* 83-17 93-99% 73-07* 
oat 97-16* 09-01% 86°81 92-49 81-89% 81:31 92:38 70-65* 
22 95-39* 07-54 §4:23 90-40 79-16 79-16 90-75 68-10* 
23 93-33 05:99 * 81-84% 88-64* 76:60 76:96 89-27 65-11* 
24 91-08 04:52 79-13* 86-43* 74:42* 74-81 87:59 62-48* 
25 88-79 02-86 76-28* —_ 84-20* 71-34 72:54 85°85 59-82* 
26 86-30 00-83* 73-46 81-89* 68°52 70:19 83-88 56-81* 
27 83-17* 79-62 65-61 67:72 81-92 53-90* 
27 88-79 898-94* 70-86 

a8 79-13* 77-16 62:67 65:11* 79-90 50-82* 

moo 84:98 96-53 67-61* 

29 74-67 98-29 74:42* 36390-11 59-57 62-48* 77-72* 47+62* 

129 81-84* 64:56* 
30 79-13* 95-39* 60-71 71-91 87-92 56-48 59-82* 75-41* 44-37* 
1 76-28% 92-91 61:44 69-19 85-65 53-19 56:81* 73-07* 40-95* 
32 73-28 90-59 57-40* 66°36 83-45 49-82 53-90* 70-65* 37-51% 
33 70-53 88-19 53-60 63-47 81-04 46-38 50:82* 68-10* 34-02* 
34 67-61* 85-85 50-12 60-47 78-48 42-85 47-62* 65-50 30-35* 
35 64-56* 83-17* 46-49 57-30 75-97 39-28 44.37% 62-48* 26-63* 
36 61-44 80-75 42-89 54-18 73-32 35-63 40-95* 59-82* 22-78* 
37 58-38 78-16 39-20 50-89 70-50 31-84 37:°51* 56:81* 18-88* 
38 55°11 * 75-47 35°47 47°51 67:65 27-89 34-02* 53-90* 14:65* 
39 51-83 72-61 31-63 44-00 64-72 23-94 30-35* 50-82* 10-57* 
40 48-44 69°77 27-70 40°57 61:74 19-95 26 -63* 47-62* 06-35* 
41 45-01 66-86 23-75 36-92 58-55 15-81 22-78* 44.37% 01-99* 
42 41-50 63°81 19-66 33-15 55-38 11-56 18-88* 40:95*  797-60* 
43 37-78 60-71 15-50 29°35 52-15 07-22 15-06* 37-51% 93-03* 
44 34:06 57-40* 11-27 25-52 48-74 02-89 10-88* 34-02* 88-45* 
45 30-22 54-11 06-93 21-51 45-28 298-36 06-91* 30-35% 83-74% 
46 26-40 50-82 02°55 17-50 41-71 93-80 02-51* 26:63* 78:92* 
47 22-47 47-42 798-05 13-26 38-05 89-09 798-21* 22:-78* 74-04* 
48 18-41 43-87 93-50 09-09 34-29 84-33 93-76* 18-88* 69-09 
49 14-30 40-29 88-92 04-72 30-49 79-49 89-24* 14-65* 64-15 
50 10-11 36°55 84-14 00-29 26°53 74:58 84.55% 10:57* 59-05 
51 05-79 32-79 79-32 295-78 22-46 69°55 79 -80* 06-35* 53-78 
52 01-39 28-86 74-37 91-21 18-45 64-43 74-96* 01-99* 48-43 
53 796-88 24-88 69-40 86-42 14-22 59-17 69:95* 797-60* 42-91 
54 92-33 20-83 64-27 81-64 09-98 53-85 64-98* 93-03* 37-55* 
55 87-64 16-69 59-09 76°75 05-46 48-52 59-83* 88-45* 31-80* 
56 82-90 12-46 53-83 71-79 01-07 43-00 54-59* 83-74% 25-92* 
57 78-06 08-18 48-53 66-72 296-58 37-32 49-34% 78:92* 20-36 
58 73-14 03-80 43-14 61-51 91-90 31-76 43 -87* 74-04* 14-32 
59 68-14 799-25 37°57 56:27 87-11 26-01 38:38 68°77 08-37 
60 63-04 94-70 32-01 50-96 82:31 20-30 32-79 63-68 02-24 
61 57-90 90-01 26-15 45-53 77-47 14-21 27:03 58:48 695-91 
62 52-63 85-27 20:56 40-10* 72:36 08-21 21-24 52:89 89-62 
63 47-29 80-38 14-73 34-48% 67-26 02-00 15-33 46-88 83-03 
64 75-48 08-75 28 -68* 62:12 195-78 09-31 39-87 76:32* 
64 41-83 44-37 
65 36-33 70-45 02-71 22-85* 56-70 89-54 03-09 37 -55* 69-12 
66 65°31 696°57 16-94% 51-32 83-11 696-77 61-16 
66 30-68 31-80* 65-40 
67 25-03 60-08 90-40 11-05* 45-97 76-62 90-33 25-92% 57-63% 
68 19-20 54-94 84-07 04-80* 40-10* 70:09 83-62 19-72 50-49 


69 13-34 49-45 HTS NIT 34-48* 63°32 LO:32% 13°55 43°55 
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Table 3 (cont.) 


P(J) QU) 
Thesis 92-41 
64:-67* 86-13* 
58-00* (ess 
SAS Oe 
44-41% 66:38* 
37-46* 59-66 
30393 52-78 
23-22% 45-86 
16-28 38°72 
09-01 31-38 
01-62 24-27 
594-08 16-48* 
86:56 08-02 
78-74 098-55 
114-78 
A 04-30 
63°36 095-48* 
55255 87-66 
47-39 79-95 
38-66 ED 
64:00 
55°79 
47-64 
39-39% 
30-98 
22-53 
14:05 
05:34 
35996 -54 
88-03* 
78:80 
69-96 
60°72 
51:49 
42-17 
32°57 
23°28 


3,0 Band 
R(J) 


28-68* 
PBK bays 


16-94* 
11-05* 
04-80* 
NOTES 
93-11 
86-13* 
19:09" 
72-92* 


66-38* 
59-14 
52-03 
44-96 


37-49 
27-73 


P(J) 
56:52 
49-86 


42:78 
Asset 
28°47 
21°28 
13-84% 
06-28 
00:05 
091-84* 


83-94 
76-13 
68°29 
60-17 
51-86 
43-94 


34-57 
23-47 


Qi) 


67-83% 
86-06 
57-63* 
76-32% 
67-83% 
60-11 
52-72 
45-50 
38-42 
31-11 
23-92% 
16-50 


08-99% 
01-75% 
01:16 
593-81 
93-20 
86-01* 


77-98 
70-06 
62-06 
53-88 
AS: Sf 
37-24 


28-89 
20-30 
ReG7s 
02-93 


493 -86* 


84-94% 


76:02 


4,0 Band 
R(J) 


07-26 


00-50 
692-61 
91-36 
83-80* 
77:20 
70-15 
63°31 
56°51 
49-26 


42-57 


Table 4. _Wave-numbers for the 4,1, 5,0 and 5,1 Bands 


424 
2,0 Band 
Q(J) R(J) 
44-14 
07°31 
38-40 
01-30 
695-14 32-89 
88:92 27-14 
82:59 21:29 
76°19 15-43 
69-64 09-35 
63-11 03-24 
56°35 697-08 
49-53 90-91 
42-75 84-29 
Geld 
35-76 
ACS He 
28-72 
64:67* 
21-50 
14-30 58-00* 
06:97 51-18* 
599-64 44-41* 
92-12 37:46* 
84-52 30-39* 
76-74 23-22% 
68-96 
60:94 
52-09 
54-20 
39-79 
45-61 
37°25 
29-05 
20:58 
12-15 
03:59 
495-00 
86-30 
77-56 
68-62 
59-61 
50:47 
41:18 
4,1 Band 
QV) 
36137-49 
35+58* 


QV) 


SE SMALAD™ 
69-98* 
68 -07* 


66-15% 
64-11% 
61-96% 


5,0 Band 
RW) 


3738609 
85°38 
84-63 
83-69 
82°85 
81:95 
80°89 
UTS} 
78:19 


76°75 
Lyon 
73-92 


PY) 


37370°89 
69:09 
O74 
65:13 
63-02 
60-81 
58°45 


55-93 
53°22 
50°70 


5,1 Band 


QV) 


36627°65 


25°87 
23502 


22:17 
20:20 
18-16 


RJ) 


PU) 
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61-38 
57-24 
52-78 
48-42 


43-94 
Sos 
34.57% 
AHS, 
24-92 


19-98 
14-95 
09°77 
04-45 
35999-34 


93-78 
88 -03* 
82-46 
76:83 
70-99 
64-95 
58:78 
52-63 
46:04 
38-89 


30°67 
48-37 
20:59 
39-17 
31-05 
23-28 
16:37* 


QU) 


Table 4 (cont.) 
5,0 Band 
Ri) P(J) 
TAO 48-02 
69-98* A723 
68-07* 42-31 
66-15* 39-15 
64-11* 36-08 
61-96* 32°88 
59-68* 29-62 
57-31* 26-22 
54:83* 22-73 
52-00 19-20 
49-32 15-36 
46°55 11-64 
43-72 07-68 
40-72 03-57 
37-58 299-56 
34-49 95-31 
31-29 90-95 
27°83 86°69 
24-10 82°51 
20-58 77°57 
16-96 72°83 
13-03 67-96 
09-00 63-13 
04-89 Sf) 
00-69 yyy 
295295 47-52 
41-96 
90-45* 36-07 
84-66 
80-24 28-03 
75-25 22-03 
70-17 15-64 
64-71 09-56 
03-05 
196-62 
56:46* 
49-27 
42-67 84:98 
36°62 76°34 
30-58 68-71* 
24-85 61-15 
19-01 54°12 
12299 47-11 
06:97* 40-02 
00-67 32-98 
194-52 25-56 
88°47 18-32 
81-98 10-80 
75:30 03.71* 
68:71* 095-88 
61:69 88-29% 
80-31% 
2230-4 


5,1 Band 
QU) 


16-20* 
14-06* 
1-562 
08-99% 
06:56* 
03 -86* 
OT Ga 


598 -30* 
D532 & 
9230" 
eelomiles 
86-01 
82-93% 
19-52 
76:03* 
72°42* 
68-82* 


65-02* 
61:04* 
Syeiibe 
52:89* 
48-80* 
44-64% 
40-24* 
35-078) 
Sigil 
26-38* 


PLES? 
L6<Si* 
test 
05-91 
00-08 


493 -86* 
86:28 


491-99 


$4:94% 
78°33 


72:19 
66-04 
59°99 
53°76 
47-77 
41-43 
34:86 
28°66 
22°18 


R(J) 


36616 -20* 


14-06* 
50% 
08 -99* 
06:56* 
03-86* 
01-16* 
598 -30* 
Deg 
92-36* 
tegen 


85-80 
82-24 
78-80 
oes 
leo! 
67-65 
63-67 
5929 


34:77 
29-99 
24:85 


L745 
10-22 


0307 3me 


497-94 
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36582-93* 
o-oo 
76-03* 
72°42* 
68 -82* 
65-02* 
61-04* 
cise 
52.00% 
48 -80* 


44-64* 
40-24* 
35-67* 
Shulah tales 
26°38* 
PHI 
16:51* 
iL eteyyfes 
06°61 

01-25 


468-90 
63-28 
56°71 


45-84 
37-53 
29-54 
23-03 
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Table 5. Wave-numbers for the 6,0, 6,1 and 7,1 Bands 


6,0 Band 6,1 Band 7,1 Band 
J QU) RV) P(J) QV) RV) PJ) Qi) RV) PJ) 
11 37857-14 
12-55-91 | 
13. 54-52 37110-30* 37586-72 
14 52-98* 08 -92* 85-31 | 
15 51-62 07-38* 83°75 
16 49-84 05-75 37114-15 81-90 
17. 48-08 04:15 13-03 80-10 
18 46-26 02:58 11-71* 37093-29 78-18 
19 44:55 00-44 10-30* 90-98 76:04 37585-99 
20 42-50 37852-98* 098-52 08:92* 88-:29* + = 73-95 84-42 | 
21. > >40-21 51-06 37830-:22* 96-44 07:38*  85-76* = 71-69 82:59 37561-54? 
22 37:95* 49-18 26-95 94-13 05:61* 83-08 69-34 80°55 58-39 
D3 seo AT ee ~ 4a 16 23-72 91-81 03-71 80-31% 66-87 78-84 55-41 
2A 32-93 45-01 20-62% 89-47 01-67 77-18 64-11 76-66* 52-59 
25. -30:22% 49.85 17-48 86:92 099-63 74-00 61°54* 74-66 49-29 
26-27-54 40:74 14-31 84-15 97-24 71-02 58-88 72-41 46-05 
Py met 71 37-95* 10-98 81-41 94-69 67°67 56-00 70-05 42-81? 
28 21-78 35-47 07:44 78-59 92-10 63:99 53-16 67°55 39-35 
29." 18°65 03:66* 75-52 89-10 60-49 50-07 65-07 35-72 
30 =: 15-33 00-20 72:36* 85-76% 56-74 47-05* 62-45 32-03 
3111 88 69-06 52-49 43-61 59-63 28-21 
32  08-52* 65°53 48:18 40-33 56-65 24-57% 
33-04-57 61-88% 36°81 53-63 20:35 
34 = 00-52 20-62% 57-74 78-13 33-19 50-54 16-44 
35 795-97 16-84 53-51 74:49 29-43 47-05* 12-21 
36 =—-91-24 12-94 777-67 48-93 70-47 35-36% 25-48 43-27 07-97 
37. —- 85-90% 43-70 21-05 38-85% 03-32 
37 08-57% 72-38% 61-21 66-45* 30-45 24-57*  — 42-81* 
3872 80°07 37-99 498-18 
38 67-61% 55-75 61-88* 25-17 19-92* 38-85* 
39-73-84 13-66 31-63% 71-64 92-67 
39 62-60* 50-32 56-92 20:25 15-31* 34-93 96-52 
40 08-57% 66-45% 
40-86-96 45-00 "51-38 14-53 10-92 31-04 91-32 
41 03:66* 64:76 61-88%  — 22:58* 
41 82-05 40-10* 08-32* 06-49 27-42 86-31 
ADs 799-10 58-27% 57-74% 16-64* 
Ve BE: 35-36* 01-35 02:14 23-58 81-40 
43 72-38% 94-61 52-39 30°80 53-30 10-77 497-58 19-49 76-21 
44 67-61* 90-36 46-19 26-07 48:93* 05-10 93°16 15-31* 71-28 
45 62-60% 85-90  — 41-18 21-32 44:51 36999-74* 88-42 11-37 66°13 
46-57-73 81-47 35-30 16-64 40-20* 94-02* 83-65 07-21 60-77 
47 —- 52-86 77-12% 29-78 11-83 36-10 88-66* 78:78 02-79 55-44 
495. 47-78 72-38% 24-02 06:95 31-63% 83-26* 73-83 498-18 50-00 
49 42-67 67-61% 18-39 02-04 27°19 77-72* 68:75 93°55 44-34 
50 = 37-51 63°20 12:71 36997-03 22:58% 72-26 63°59 89-22 38°76 
51 32+ 17 58:27* 06-93 91-78 18-05 66°57 58-41 84:25 32-98. 
52 26-88 53-58 01-15 86-70 13-19 60:96 52-98 79-49 27-10 
53-21-32 694-94 81:25 08:32* 54-88 47-57 74-53 21-26 
54. 15-75 89-05 75-94 03-35 49-04 41-97 69-41 15-11 
55 72:91 36997-70 43-02 
55 09-82 70-17 99-74% 36-50 64:10 08-95 
56 = 05-58 65-96 36-66 
56 33°36 62-75 94-18* 30°63 59-04 02°83 
57 699-07 59-62 29-70* 
57 28-05 54-58 88:66¥ 31-99 24-87 53°64 396-39 
58 92-78 22-51 64:16 53-72 83-23% 24-99 18-75 47-95 89-96 
59 ~—- 86-90 16-99 57-64 47°86 78-01%  — 18-53 12-66 
60 80-72 11-32 50-71 41-85 72-68 11-91 06-51 
61 74-28 44-13 35-74 66-89 05°39 00-48 
62 68-10 37-54 29:70 61-24 393-83 
63 61-83 23-37 55-44 


64 BSED 17-14 49-71 
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6,0 Band 
R(J) 


QY) 


48-51 
41-76 
35-07 
27-97 
21-03 


13-87 
06-50 
599-23 
OT-79 
84-42* 
76-66* 
68-66 
60-75 
52°59% 
44-61 


36°36 
2 -A2* 
18-24 


PJ) 


Table 5 (cont.) 


6,1 Band 
QV) RV) PU) 
~ 10°73 43-65 
37-56 
31-63 


427 


7,1 Band 


QU) RV) PY) 


§4. DETERMINATION OF ROTATIONAL CONSTANTS 
The rotational constants for the ground state have already been determined 


(Barrow 1946). 


The values of B” and D” were found to follow the equations 


B,” =0-3036, — 0-0014,(0" +4); D,” =[2-01, —0-01(v" +4)] x10-. 


In order to calculate the rotational constants for the upper perturbed levels 


we have used the method introduced by Geré (1935). The following expressions 
are formed: 
+ QU-1)-QW) 
2S 
R(J —2)—RV —1)+ PY) -— PV +1) 
ee 
If the left-hand sides of these expressions are plotted against J?, AB is obtained 


as the intercept at J =0. Only the unperturbed regions, where the points lie on 
almost horizontal lines, are used for the extrapolation. These values of AB and 


Be Bie IND 


Pie BOD ID eee (2) 


Table 6. Rotational Constants for the p1II State 

v Bobs Beate Dops Deaic 

0 0:2653, 0:2653, 2°85 x1077 2:90 x 10-7 
1 0:2632, 0:2632, 2:90 2°87 

2 0-2610, 0-2610, 2°87 2°84 

3 02589, 02589, 2-82 2-81 

4 0-2567, 0:2567; 2°73 2:78 

5 0:2546 0:2545, 2°75 

6 0-2524 0:2524, 2°72 

7 0-2503 0-2502, 2:69 


of AD were used to obtain values of B’ and of D’ with the aid of values of 
B’ and D” already known. The perturbations in v’=5, 6 and 7 are, however, so 
considerable that no reliable values of D’ can be obtained for these levels. 
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The constants are given in table 6. The calculated values are from the 

smoothed expressions 
B,| =0-2664, —0:0021,(v' +4), D,’ =[2-91 —0-03(v' + $)] x 107. 

We estimate the errors in the B’-values to be about +0-00005 for the levels 
v' =0 to v' =4: for B,’, By’ and B,’ they are about +0-0001 cm~. 

With the values for w,” and w,’ in §5 we obtain, from Kratzer’s relation 
(D =4B3/w2), D’=1-99 x 10-?, D’ =2:87 x 10-7, in satisfactory agreement with 
the experimental values. 


§5. DETERMINATION OF VIBRATIONAL CONSTANTS 

The origins of the bands analysed have been given in table 1. The vibrational 
-constants have been calculated in the usual way, omitting the perturbed origin 
of the 7,1 band from consideration. We obtain w,” =749-69, w,’ =513-12, 
x,"w,” =2:58, x,/w,’ =2:93, v,.=35026-86. With these values the calculated 
origin of the 7,1 band is found to be 37591-72; this is 4-69 cm™ lower than the 
value observed. 

§6.. PERTURBATIONS 

It is found that all the upper vibrational levels are perturbed: the lower 
levels are unperturbed. The perturbations have been detected by forming 
the expressions (1) and (2) and plotting them against J (fig. 2). Information 
about the perturbations is summarized in table 7 and in fig. 3. ‘The figure shows 
the upper vibrational levels in a plot against J(J+1): circles, squares and 
triangles indicate the positions where the perturbations culminate. 


The Perturbing States 

The appearance of the perturbations shows that they are caused by 
interactions of two perturbing states, a triplet one and a singlet one. These 
states are here called e and f. Both states lie lower than the p1II state. It is not 
possible to decide whether the state f is 'II or 1A. It cannot be 12 as both the 
A-type terms of D'II are perturbed. State eis probably ® , asin the CO spectrum 
(Coster and Brons 1934, Budd and Kovacs 1938). There are three interactions 
ata *®2 .1II perturbation: the c terms in II are perturbed by the °X; component 
and the d terms by the °X;_, and *%;,, components. This means that at a 
3> .1II interaction one perturbation may be found in the Q branch at a certain 
J value, two in the R and P branches, one at a lower, the other at a higher J value 
than that for the Q branch. In fig. 3 the *X; .‘II perturbations are indicated by 
squares, the *2;,,.'Il perturbations by circles and the f.1II interactions by 

triangles. 

Determination of Constants for Perturbing States 


At several perturbations extra lines have been found. If extra lines are 
found at least at two consecutive J values, it is possible to determine the 
B and C values of the perturbing states (Kovacs 1937, 1951: for practical 
examples see also Lagerqvist, Lind and Barrow 1950 and Hultin and Lagerqvist 
1951). ‘The ordinate of an intersection point in fig. 2 gives, for a singlet-singlet 
perturbation, 3(B,+B’)—B”, where B, is the B value of the perturbing state. 
Kovacs (1951) has shown that if the terms of a 3X state can be written as 


Fyy,=C+Bix(J+1)(J+2), Fy =C+BisJ(J +1), 
alee =C+ Bss(J — 1)J 
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Fig. 2. [Q(J—1)—Q(J)]/2J and [R(J—2)—R(J—1)+P(J)—P(U +1)1/4J plotted against J. 
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then the middle component of a *X.1II perturbation may be treated in the same 


way as for a singlet-singlet interaction. The B values for the other components — 


are 


B,s(P) = {J/(J + 1)}B,(P) for K=J +1, 

Bis(P) = {J /(J — 1)}B,(P) for K=J —-1. 

The values of C (the separation at J =0 between the perturbed and perturbing 
states) are given by 


C(Q) = —2%) + B’- LJ -1)[QV) + QJ) -(F + DIQU —1) + O'Y — 1} 


op sonnet (4) 
CP) = —219 + 2B" + B’—-i{(J —1)[/PV +1)+ P'(J+1)4+RJ-1)+R'U-D] 
—(J +1)[PV)+P'(J)+ RJ —2)+R'(J—-2)]}. we. (5) 
An approximate value for C can be found from 
C{Q)=BiUY 1) ba ee (6) 
where J* is the J value where the curves in fig. 2 cut each other. 
40 
«10° 
39 
V=7 
38 
3 6 
E 
> 5 
37 
4 
3 Zs 
36 be 
p= 2 = ‘7 ‘TT 
e y+, 170 
| a f'Tl 
35 0 
0 l 2 3 4 mar 6 7 8 9 io if 7 


Fig. 3. Perturbations in the state p 4IT of SiS; term values plotted against J(J +1). 


The results of the treatment of the perturbations are collected in table 7 
The B, values given there are the means of the observed values (the number éf 
observed values is given in brackets). Only at one °X” perturbation has an 
intersection heen found in the R and P branches, namely at /~40 in the 6,1 band: 
the values of C+ in table 7 are therefore those calculated from the Q fines: 
using eqns. (4) and (6). It has not been possible to determine the absolute 
numbering of the vibrational levels of the perturbing states. The lowest observed 
levels have therefore been given the vibrational quantum numbers e and f. 
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The perturbation between eX and p!II in v' =0 at J~66-5 in the Q branch is 
very weak, but it is clear enough to be detected. This can be seen from the 
behaviour of the Q branch of the 0,0 band in fig. 2. The perturbation can also 
be traced in the bands analysed earlier, as may be seen from the plot of the mean 
value of [Q(J —1)— Q(J)]/2J for the 0,1, 0,2 0,3 and 0,4 bands given in fig. 2. 
The corresponding perturbations in the R and P branches have not been 
detected: presumably the perturbing matrix elements are too small. We have 


Table 7. Summary of the Perturbations 


v’ J* By Us Cst vo Perturbing 
RE Q RP eqn. (4) eqn. (6) state 

0 — 66-5 — i e€ — SIMs EDS 

de 45 bie? 54 O20 11 e-- 1 Sb Sote 30040/ nme ie 

1 — +>100 — a e+2 — — DR 

PO NS, 27°8 30 0:2092(2) ar 2a S09 5/3 5950-4 en 

2 — 03 — — e+3 — 36364-3 am 

Seedy. 83°5 (86-5) — e+4 — 36769:0 9 3 

B05 °5 70°8 74 0:2094(2)'— e==5 37144:9 371553 25 

4 — 82:5 _ _ i — oS 1TT or 1A 

S51 56°5 aus 0:2056(4) e+6 37546-1 S7545:0 eee 

pee /() 70 — — f+1 _- — HIE exe Gah 

oR Bos S/S 40 0-2013(3) e+7 37940-1 3/1932: a ie 

Gu 20-5 56:5 — 0-224(2) f+2 = — TA eye YN 

7 
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also been unable to find the 2X" .1II perturbation between the levels v' = 1(p1I1) 
and v,=e+2/(®2 ), which also seems to be too weak for detection. It should 
be situated at /~102. 

The B values for e?2” may be summarized: B,,,=0°213,—0-001,v. For 
the state f only one value has been computed: B,,.=0-22;. Since the absolute 
v-numbering is unknown, the vibrational constants must be left as follows 


e°2: -o, = 415, 'x.w, = 179, =f 7 w~450 em. 


§7. CONCLUSION: COMPARISON OF SiS WITH CO 

The information at present available about the electronic states of SiS is 
summarized in table 8, and in table 9 a comparison is made between SiS and CO, 
the first member of this group of molecules. It is clear that the analogy between 
the states of these molecules has been considerably strengthened by the results 
of the present work. ‘The two states, e*& , both of which are known only through 
their interactions with other states, are certainly very similar in properties. 
Much less is known about the perturbing states f!II or 1A, both in CO and in SiS, 
but from what is known they too appear to be similar. ‘Thus four of the six 
known states of SiS appear to have their counterparts in the term-scheme 
for CO. It seems however that it may be dangerous to attempt to take the 
analogy much further at the moment. ‘T'hus, a third perturbing state, I'D’, is 
known from its interactions with a'll in CO, but no corresponding state has 
yet been found in SiS. It is of course possible that such a state will be revealed 
by perturbations in the D—x system of SiS lying outside the range of the present 
measurements, although such work is likely to be very difficult on account of 
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the weakness of the bands concerned. More striking is the difficulty of finding 
a counterpart to the E state of SiS in CO. The behaviour of the force constants 
has led to the suggestion that this state may correspond either to d*II, or to a 
a’8+ of CO (Vago and Barrow 1946). Neither identification can, however, 
be accepted very readily, for the E-x system seems to be stronger in absorption 
than would be expected for a forbidden transition. It may therefore be that a 


Table 8. Summary of Constants for SiS 


State Ve We XeWe Be ie DOB 102s 7. (A) D(ev) 
E 41923:5. 403:5, 1:40*<0-30 — — ——. S193 12a 
f4IT or 1A <37150 (450) — 20-224 <2:24 — 
ele <35310 3415 1:9 20-214 0-001; _— — <2:30 — 
pil 35026:8,1513-15 92°95 0:2664, 0:0021, 2-91 2 2 058om2 tei 
at (15000) (650) — — — — 
xi 0 749:6, 2:58  0-3036;0-0014, 2-01; 1 1:928, 6:4,, 


* —0-032,(v-+4)?—1-1,; x 10-°(v-++ $)° (Thomas 1947). 

+ By extrapolation of the vibrational levels in state E, and assuming that states x, D and E£ all 
dissociate into ground-state atoms, Si(?P)+S(?P) : the multiplet separations in these atomic states 
introduce a maximum uncertainty in the values of Dy of about 0:1 ev. 

t This state is very uncertain at present. It is included here on the assumption that the visible 
system attributed to SiS arises from a transition E—a. No satisfactory analysis of this system has 
yet been made, however, and the visible bands may well arise from a transition between two states 
of SiS, neither of which is yet known to take part in other transitions. 


Table 9. States of SiS compared with the Lowest States of CO 


ere) SiS 
State tig wax B/Bx State De w/wx B/Bx 
d *; 61785 0-52 0-65 
e 8 <64800 +2049 +0-64 efn- <35140) <30:55 i) 20-70 
a’ 3+ (55380) 0:56 0-69 
a *II, 48474 0:80 0:87 
E? 41751 0-54 eat 
B1Dt 86918 0:96 1-02 
TiS <66380 30-48 >0-75 
MTT or 1A < 66230 = < 0:83 f1II or 1A <37000 3060 >0-74 
A2II 64747 0-70 0-83 p4UI 34909 0-68 0-88 
(a) ? (15000) (0-87) =* 
on 0 1-00 1:00 x yt 0 1-00 1-00 


The data for CO have been taken from Herzberg (1950), except for the constants for the 
perturbing states. The figures for e*}~ and 1'D~ are derived from Schmid and Gerd (1937). 
The existence of the ‘IT or 1A state of CO at v<66230 is based on the observation of a perturbation 
near the beginning of v=1 in a'I] by Geré (1935) and Schmid and Gerd (1935). 


real difference between the electronic schemes of the two molecules is appearing 
here, and that the counterpart to the E state of SiS is to be found, not in CO 
but, for example, in the isoelectronic molecule P,, which has an excited 1E+ state 
of weaker binding than in the ground state. A rotational analysis of the 
E-x system of SiS which has been begun by one of us (G.N.) should throw 
light on this question. ‘Too little is known about the only state of SiS remaining, 
the lower state, a, of the visible bands, to justify much consideration now. The 


most likely interpretation, on the basis of table 9, is that it may correspond to 
a°lPeot CO. 
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Dissociation of Molecular Hydrogen Ions in the Cyclotron 


By Ky Ey AD EFFAT 
Department of Physics, The University of Birmingham 


Communicated by F. H. Fremlin; MS. received 28th Fanuary 1952 


ABSTRACT. The dissociation of H,* ions by collision with gas molecules greatly reduces 
the current of these ions that can be obtained from a cyclotron. ‘The current has been 
measured at two distances from the ion source corresponding to energies of 9 and 18 Mev 
and with various pressures of air and argon in the tank; the loss of ions in expanding from 
the smaller to the larger radius is somewhat less than is predicted by Salpeter’s calculations 
of the dissociation cross section. Assuming this cross section to be inversely proportional 
to the kinetic energy of the ion, the product cE is found experimentally to be 
(1-9+ 0-4) x 10-16 cm? Mev per ‘atom’ of air and (7:0+1-4)x10-'6 cm? Mev per atom 
of argon. 


HE Birmingham 60-inch cyclotron (Fremlin and Gooden 1950) is a fixed 

machine which is normally used to accelerate molecular hydrogen ions H,*, 

deuterons D+ or alpha-particles He?*, using nearly the same magnetic 
field in each case. ‘The beam of 20 Mev molecular hydrogen ions is equivalent 
when incident on a target or when emerging through a window to a beam of 
10 Mev protons. The acceleration of H,* ions has therefore the advantage that 
the proton beam is obtained with the same magnetic field as for deuterons and 
alpha-particles ; the re-shimming of the magnetic field which would be necessary 
if protons were to be accelerated directly is theretore avoided. 

With H,* ions the circulating beam current at full energy is, however, some 
10 to 15 times less than would be expected from comparable observations with 
deuterons, and this suggests, as pointed out to the author by Dr. J. H. Fremlin, 
that dissociation of the H,* ions occurs as a result of collisions with gas 
molecules during acceleration. 

Since the energy required to dissociate a molecular hydrogen ion is only a 
few electron volts, the dissociation cross section cannot be studied by the usual 
method of interposing foils into the beam. In collisions between high-energy 
molecular ions and gas molecules in which such small amounts of energy are 
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transferred, a proton and a neutral hydrogen atom or two protons each of half 


the incident energy will be formed. The protons will then move in circular _ 


orbits of radius half that of the orbit of the incident H,+ ion as shown in fig. 1. 
This should cause most of the protons so formed to strike the ion source situated 
at the centre of the cyclotron. Visual observations through the optical system 
arranged so that the inside of the vacuum chamber can be viewed from the 
control room show that the ion source tip becomes white hot when H,* ions 
are accelerated, but not during the acceleration of either deuterons or 
alpha-particles. If the gas pressure in the cyclotron tank is increased, more 
H,+ ions are dissociated and the ion current therefore decreases. For ions of 
energy above a few Mev dissociation is much the most probable cause of removal 
of ions from the beam owing to the small energy transfer required; and the 


A 
“Proton from Dissociated 
H,* lon at A 


Neutral 


Fig. 1. Diagrammatic representation of dissociation of H,* ions in the cyclotron. 


dependence of ion current on gas pressure can therefore be used to measure the 
dissociation cross section. A calculation made by Salpeter (1950) shows that 
the cross section for dissociation of an H,+10n of energy £ is expected to have the 
form o=o,/E. 

If then the ion current is found under steady conditions at a radius R, from 
the centre of the cyclotron and then at a greater radius R,, the two currents 
I, and J, will be given by 


In 


Hf -6,(R,—R 

1 _ x Colhe  », ali (1) 
Te AE 
where p is the gas pressure in the tank in microns, and J, and /, are the ion currents 
at the radii R,, R, respectively, o (cm? Mey) is a constant for a given gas which 
equals numerically the dissociation cross section for an H,* ion at an energy 
of 1mev, and AF 1s the energy gain per revolution in Mev, assumed to be 
equivalent to an acceleration of twice the peak voltage between the dees; 
K =4-45 x 10" for a monatomic gas. 
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lon currents were measured calorimetrically so that uncertainties due to 
radio-frequency pick-up and secondary electron emission were eliminated. 
Measurements were carried out at radii of 18 in. and 25! in. from the centre 
of the cyclotron, corresponding to nominal H,+ ion energies of 9 and 18 Mev 
respectively. Over this region the loss of ions is mostly due to dissociation 
other processes, such as electron capture, scattering, etc., are not important 
because momentum transfer in ionizing collisions is small and nuclear collisions. 
are rare. 

Dry air or 99%, argon was admitted to the cyclotron tank through a needle 
valve to vary the pressure between the ultimate vacuum, 3-4 x 10->mm, 
with the cyclotron running and 10-*mm Hg, which was the highest pressure at 
which the accelerating voltage could be maintained on the dees. The pressure 
was measured by a B.A.R. ionization gauge connected by a wide tube to the 
cyclotron tank and located in a position where the effect of the pressure gradient 
due to the pumping manifold and baffling is not important. Pressure changes 


H,* lon Current Ratio /,//, 


Ratio with 
Deuterons 


lon Gauge Plate Current “A (pressure) 
0-854 > 10°mm for Air 


Fig. 2. Plots of the H,* ion current ratios at 9 and 18 Mev against tank pressure. 


of 10-*mm could be obtained with the leak system and detected by the 
microammeter (sensitivity of 10-7 amp/div.) measuring the positive ion current 
of the ion gauge. Absolute pressures were obtained using conversion factors 
of the ionization gauge (Riddiford 1951), The r.m.s. voltage between each 
dee and earth was measured by diode rectifiers connected by quarter-wave 
transmission lines to circular probes placed near the dee surfaces. 

The results of typical runs for air and for argon are shown in fig. 2, and 
confirm the linear relation between In /,//, and p predicted by eqn. (1). The 
intercepts obtained by extrapolation indicate that, in addition to the losses due 
to dissociation, there are losses which are independent of the pressure. ‘These 
losses appear to be a function of the operating conditions of the machine, and 
may be caused by deceleration of the ions by excessive phase change since it is 
found that the ratio of the currents obtained at the same radii with deuterons 
agrees reasonably well with the intercept obtained by extrapolation of the curve 
for air. In the case of argon, the background pressure is due to the gases present 
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in the tank when the argon is introduced; extrapolation is therefore carried out 
with a slope corresponding to the air cross section as indicated by the dotted 
line on the curve for argon. This also gives an intercept in agreement with the 
ratio of currents obtained with deuterons. The absolute average cross sections 
deduced from independent sets of observations and the corresponding values 
of the mean free path for dissociation at a pressure of one micron are compared 
with the two estimates given by the theory (Salpeter 1950) in the table. 


Ca 9 (cm? Mev) (Mean free path)/E 
per atom x 1018 (cm/Mev) for p=1 micron 
Theoretical Experimental Theoretical Experimental 
Air 220-350 190 + 20% 41-66* 75 2073 
Argon 990-1300 700 + 20% 22-31 40+ 20% 


* In the preparation of the table given by Salpeter (1950) some factors of two were inadvertently 
omitted. The results quoted in Salpeter’s letter for (mean free path)/E for hydrogen, nitrogen 


and oxygen should be divided by 2. 


The theoretical values for air is obtained from the values given for nitrogen 
and oxygen. 

The difference between the cross sections is about 30° with an estimated 
experimental error of 20°, which is reasonable in view of the uncertainties of 
the theory. The ratio of the cross sections of air and argon agrees with the 
theoretical value to within 5%. Previous investigations have only been carried 
out at low energies, up to a few kev, and were not intended for measurement of 
the dissociation cross section. Estimates of 5 to 300 collisions per cm path, 
at 1mm pressure for 1160 ev H,* ion were given by Lamar et al. (1935). These 
values cannot be compared with the present results because the inverse 
proportionality between o and F will not be valid in this energy region. 

It is clear from these experiments that good vacuum conditions are of great 
importance for the acceleration of molecular ions. In this cyclotron, with a 
pressure of 4x 10->mm, 70% of the H,+ ions reaching an energy of 1 Mev are 
later lost by dissociation. ‘This loss is reduced to 3% at a pressure of 10-6mm, 
which can be obtained in large vacuum systems. It is interesting to note that 
in the case of a synchro-cyclotron, the vacuum required for efficient acceleration 
of H,* ions would be better than is obtainable with present techniques. In a 
typical 200 Mev synchro-cyclotron, with a total path length of about 100km 
50% of the H,* ions would be lost by dissociation at a pressure of about 
2x10-?mm. Ina1000 mev proton synchrotron where the path length is much 
greater the loss of particles would be virtually 100%. 
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ABSTRACT. An experiment is described in which the momentum and specific ionization 
of cosmic-ray mesons were simultaneously measured, the momentum by means of a cloud. 
chamber in a magnetic field, and the specific ionization by rectangular proportional counters. 
The difficulties inherent in the use of proportional counters for such measurements are 


discussed. 
‘The experimental results are in satisfactory agreement with theoretical predictions of an 


increase of specific ionization for very energetic -mesons. 


§1. INTRODUCTION 
1.1. Theoretical Background 

HE average specific ionization produced by fast charged particles in their 
| passage through matter has been discussed by Bohr (1915), Bethe (1933) 
and Bloch (1933 a, b). They conclude that the ionization should decrease 
with increasing particle momentum until a minimum value is reached at a 
velocity given by v/c=$=0-96; this minimum should be the same for all 
particles with equal charges. When the momentum is raised still further, the 
ionization should slowly increase. For convenience this effect will be referred 
to as the ‘logarithmic increase of ionization’. It is a relativistic effect, due to 
the Lorentz contraction of. the electric field, and the consequent increase in 

the maximum distance at which a particle is able to ionize. 

Swann (1938) and Fermi (1940) pointed out that polarization effects in the 
medium should screen the more distant atoms from the electric field of the 
particle. As a result, the ionization should not increase indefinitely as 
B=v/ce+1, but should reach a limiting value, and then remain constant. 
Fermi’s one-electron theory has been generalized by Halpern and Hall (1948). 
For low energies their results agree with those of Bethe and Bloch, but 
differences between the two theories become important when the particle 
velocity exceeds the velocity of light in the medium. ‘The differences are, 
however, negligible when .-mesons of momentum less than about 10'°ev/c pass 
through the proportional counters used in the experiment described below. 

The polarization effect (which is more usually referred to as the ‘density 
effect’) should be particularly important in solids and liquids. Detectors such 
as crystal counters and photographic plates are therefore rather less suitable 
for investigations of the logarithmic increase than cloud chambers or gas-filled 
counters. 

So far we have referred to the average ionization, which is the quantity 
usually considered theoretically. There are good reasons, however, for 
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preferring to deal in experimental work with the most probable ionization. 
The statistics of the ionization process were first considered by Bohr (1915) and 
Williams (1929). They pointed out that large fluctuations in the rate of loss of 
energy must occur from point to point along the track of a particle. These 
fluctuations arise partly from the Poisson distribution of the number of collisions 
per unit length, and partly from the widely different energies which can be given 
to the secondary (8) electrons at each collision. The probability of any particular 
energy-loss occurring in a given track-length has been calculated by Landau 
(1944). His results show a broad distribution for which the most probable 
ionization (the peak of the distribution) shows the same kind of dependence on 
v/e as the mean ionization*. Either value could equally well be used to verify 
the existence of the logarithmic increase. From the experimental point of view, 
however, the determination of the mean presents several difficulties, of which 
saturation of the detector or amplifier by large pulses is an obviousexample. ‘These 
difficulties are avoided by measuring the most probable specific ionization. 


1.2. Previous Experimental Results 

A considerable number of attempts have been made to observe the 
logarithmic increase and the density effect experimentally, with rather 
contradictory results. In experiments using cloud chambers the ionization is 
measured by drop-counting. Occasional high-energy collisions cause dense 
clusters of drops which cannot be counted accurately, and have therefore to 
be excluded from the results. However, there is no serious objection to this 
procedure provided the most probable, and not the average ionization is being 
measured. Corson and Brode (1938), Sen Gupta (1940) and Hazen (1945) 
made such measurements on electrons, and all detected an increase of ionization 
with momentum. Similar experiments on p-mesons by Sen Gupta (1943) 
and Hazen (1945) suggested that the rise was smaller than had been predicted. 
On the other hand, Ghosh, Jones and Wilson (1952) obtained results which 
agreed very well with theory. ‘They were able to extend their measurements to 
sufficiently high momenta (3 x10ev/c) for the polarization effect to be 
detectable. 

Low-efficiency Geiger counters, the efficiency of which depends on the 
primary ionization, have been used by Hereford (1948) to measure the 
ionization of electrons. He found a small rise as momentum increased. 
Experiments using proportional counters to measure the ionization of mesons 
of momenta up to 5 x 10!ev/c have been made by Goodman et al. (1951). They 
failed to detect any increase above the minimum, and concluded that mesons 
behaved in this respect differently from electrons. 

The results obtained with gaseous detectors therefore seem to be in favour 
of a logarithmic increase for electrons, but leave room for doubt in the case of 
o-mesons. 

The evidence provided by solid detectors is also not unanimous. Whittemore 
and Street (1949) measured the energy-loss of mesons in a crystal ionization 
chamber, and found a small increase of ionization with momentum. It was less 
than expected, however, even after allowance had been made for the density 
effect. Bowen and Roser (1951), using an anthracene scintillator, failed to 


* It is now known that Landau’s theory underestimates the width of the distribution for heavy 


gases. Blunck and Leisegang (1950) show that, by taking the shell structure of the atom into account, 
a broader distribution can be obtained. 
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detect any change in ionization for j.-mesons between 3 x 108 and 3 x 10°%ev/c. 
They concluded that the density effect prevented any rise from taking place, 
although their results were equally compatible with the view that the logarithmic 
increase does not occur for p-mesons. Pickup and Voyvodic (1950) used 
photographic plates to measure meson energies by scattering and ionization by 
grain counting. They found a 10% rise in ionization above the minimum, and 
evidence for a density-effect ‘plateau’ at high energies. 

In the present experiment we have used a proportional counter to measure 
the energy-loss of cosmic-ray mesons, and a cloud chamber to determine their 
momentum. Our results show an increase which is in reasonable agreement with 
the theoretical predictions. The work was carried out at the Aiguille du Midi 
laboratory (3 650 metres) during the course of some experiments designed to 
measure the vertical slow proton intensity. ‘These experiments will be described 
elsewhere. 


§2, APPARATUS 
2.1. General Description 


The experimental arrangement is shown in fig. 1. A simple telescope 
defined the cosmic-ray beam, which passed through a double proportional 
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Fig. 1. (a) Section perpendicular to Fig. 2. Block diagram of circuit used for 
counter axes. (5) Section parallel to recording proportional counter pulses. 


counter axes. 


counter and then through a cloud chamber. A filter consisting of 15cm of 
lead was placed above tray A, in order to cut out the soft component. Below 
tray C was a fourth counter tray D, the two being separated by 8 cm of lead. 
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The cloud chamber was expanded and the proportional counter pulses — | 


recorded whenever a triple coincidence (A+B+C) occurred. In addition 


a telephone register indicated when events of the type (A+B+C—D) took _ 


place. This information was used to estimate the number of particles stopped 
by the absorber. 

The cloud chamber, of dimensions 56cm x 16cm x 15-5cm, was situated 
in a very uniform magnetic field of about 2500 gauss. This field, which was 
produced by a large air-cored electromagnet, enabled momenta of up to 
5 x 10%ev/c to be measured. The cloud chamber, and the associated apparatus 
(not shown in the figure), were designed by Nageotte, and are described in a 
paper by Bastin et al. (1950). 


2.2. Proportional Counters 

(a) Construction. ‘The double counter used throughout the experiment 
consisted of two rectangular chambers, placed one above the other in the same 
envelope. The common dividing wall was made of thin (35 inch) brass, 
perforated to allow the same gas filling to be used for both counters. The gas 
used was a mixture of 40-6cm Hg of krypton and 4:-2cm Hg of methane 
(at 18°c). The depth of the counters in the direction of the cosmic-ray beam 
was 3-19cm in one case and 3-42 cm in the other; the internal width was 7-2cm 
and the useful length 25cm. The. wires, placed centrally’'in the two 
compartments, were of 0-004 inch diameter tungsten. The counters were 
operated at 2200 volts. 

The rectangular shape of the counter, chosen to give equal path-lengths for 
all particles, led to some difficulties in the preparation of the filling gas. It was 
found necessary to work with very pure gases in order to avoid recombination 
in the low-field regions near the corners. 

To test that a filling was suitable, calibration x-rays (Rothwell and West 
1950) were introduced into each counter through five thin windows. ‘Two of 
the windows were close to a low-field region, while the other three were on the 
centre-line of the counter, where the electric field was strong. A good filling 
gave substantially identical calibrations for each of the five windows. Impurities 
caused the two ‘low-field’ windows to give x-ray peaks that were too low and 
too broad. 

With argon practically identical curves were easily obtained for all five 
windows, but-with krypton more difficulty was found. The filling chosen for 
the experiment, using krypton that had been purified by standing over heated 
calcium for 24 hours, gave a 4%, difference in most probable pulse-size between 
windows, and a 9%, half-width at half-height for the 8 kev K x-ray of zinc. (A good 
cylindrical counter gives a half-width of 74%.) ‘This was considered admissible 
in view of the large fluctuations known to exist in the ionization of cosmic-ray 
particles, but, as a precaution, the counter telescope was arranged to include 
only the central 5 cm of the proportional counter width. 

(6) Routine calibration. X-ray calibrations were carried out regularly 
throughout the course of the experiment, using a series of artificially prepared 
radioactive sources. Lines at 8-1, 10-5, 20-2 and 41-6kev, given by zinc, 


selenium, palladium and gadolinium respectively enabled non-linearities in the, 


counter and amplifiers to be corrected, and the output pulse-height converted 
to an energy-loss. ‘T’his set of calibrations enabled us to measure energy-losses 
of up to eight times the most probable value for minimally ionizing mesons. 
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‘The mean value of the calibrations for the five windows was used, in order to 
reduce errors due to variations of sensitivity within the counter. 

(c) Recording of counter pulses. Figure 2 shows the block diagram of the 
circuit used to record the proportional counter pulses. ‘The signals from the 
head-amplifiers were separated in time by a delay-line operating on one pulse 
only, then mixed and amplified. The amplified pulse was fed through a further 
delay-line to the Y plates of a single-sweep recording oscilloscope, which was 
triggered and photographed whenever a coincidence pulse was received from 
the counter telescope. 


§3. EXPERIMENTAL ERRORS 
3.1. Errors associated with Ionization Measurements 


There were three sources of instrumental error that were not eliminated by 
calibration with x-rays: 

(i) Statistical fluctuations in the counter multiplication processes gave rise 
to a random error with a standard deviation of about 84%. ‘This figure was 
calculated for a release in the counter of 4 kev (roughly the most probable value 
of the ionization minimum) by assuming that the uncertainty would be the 
same as that found for an x-ray calibration line of the same energy. 

(ii) ‘The uncertainties in calibration, due to variations of sensitivity within 
the useful volume of the counter, gave a random error of standard deviation not 
greater than 3%. 

(iii) The r.m.s. error in reading the photographic record of the oscilloscope 
trace was estimated to be not more than 5°, by comparison of the results obtained 
by two independent observers. 

Since these errors were independent, the overall r.m.s. error of a single 
observation was taken as the root of the sum of the squares, that is, approxi- 
mately 10°. It will be shown later that random instrumental errors of this 
magnitude have only a slight effect on the accuracy of the results. 


3.2. Errors in Momentum Measurement 


The errors in momentum measurement vary with the magnetic field and 
with the quality of the cloud chamber photograph. ‘The principal error, 
however, arises from turbulence, which causes an apparent curvature of the 
tracks even in the absence of the magnetic field. ‘This error, and any errors 
due to multiple scattering, were estimated by taking about ten per cent of all the 
photographs without field, and measuring the curvature of the tracks. It may 
be concluded that, for the majority of the observations, the uncertainty in the 
momentum p is such that 1/p is known to better than +0-2, if p is in units of 
10%ev/c. The highest value of momentum that could be regarded as measurable 
was therefore taken as 5 x 10% ev/c, and all particles of momentum greater than 
this were placed in the ‘unresolved group’. At 10®ev/c the maximum error 
was +25% and —17%; at 5x 108 ev/c it was +11% and —9%. 

Due to the asymmetry of the momentum error the momentum will, on the 
average, be overestimated. If the logarithmic increase exists this error will 
cause the most probable ionization predicted on the basis of the momentum 
measurements to exceed the observed value by an amount which will be 
negligible below 2 x 10% ev/c, and less than 6% at 4x 10®ev/c. ‘The momentum 
error will therefore cause the observed rise to be a little less than that expected 


theoretically. 
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$4. TREATMENT OF RESULTS 
4.1. Rejection of Unwanted Particles 


Almost all the particles recorded by the apparatus were either -mesons, 
m-mesons or protons. As we shall show later, the 7-mesons comprised only a 
small fraction of the total beam. Since their predicted ionization curve is very 
close to that for j.-mesons, their presence is unimportant. On the other hand, 
the protons, which differ appreciably from mesons in their ionization properties, 
were sufficiently numerous at the altitude of the laboratory to make their exclusion 
essential. 

The simplest way of ensuring this was to ignore all positive particles 
possessing more than the momentum at which protons could just traverse the . 
telescope. This was approximately 5 x 10%ev/c. Below this momentum there 
was no objection to using both positive and negative particles, since the 
ionization of mesons is independent of their sign. 

At momenta greater than about 5 x 10®ev/c the tracks were so little curved 
that it was impossible to decide whether the particles were positive or negative. 
There was therefore no way of selecting only mesons from the ‘unresolved’ 
group’ of very energetic particles. 


4.2. Rejection of Multiple Events 


Approximately two-thirds of the cloud chamber photographs showed single 
particle tracks, the remainder being confused by showers or the accidental 
passage of two particles during the sensitive time of the chamber. Proportional 
counter observations were accepted only if the particles were definitely single. 
It is realized that this method of rejecting multiple events is not completely 
certain. For instance, two diverging particles could pass through the 
proportional counters and give rise to only one track in the cloud chamber. 
It is difficult, however, to imagine a mechanism whereby this kind of event 
could imitate a logarithmic increase. For this to happen the angular divergence 
of the associated particles would have to be greater the greater their energy. 
In fact the opposite is known to be true. ‘The rise we have observed cannot, 
therefore, be explained in this manner. 


4.3. Fluctuations in the Ionization Process and their Effect on the Treatment 
of Results 


The existence of fluctuations in the rate of ionization of fast charged 
particles has already been briefly mentioned. ‘The probability distribution for 
the ionization liberated in a layer of matter having a thickness much smaller 
than the range of the particle has been calculated (with the aid of some 
simplifying assumptions) by Landau (1944). He finds that on the low-energy 
side of the peak the curve should fall rapidly to zero, while on the high-energy 
side it should fall much less sharply and show a pronounced ‘ high-energy tail’. 
The values predicted by his theory for the ratio of the full width of the 
distribution (at half-height) to the most probable ionization range from 0-25 to 
0-45 for many commonly used counter gases. However, it has since been pointed 
out by Blunck and Leisegang (1950) that for heavy gases the widths should be 
even greater than those predicted by Landau. This is confirmed by some 


Ionization by Energetic Cosmic-Ray Particles 443 


recent work (Rothwell 1951) on the ionization losses of fast electrons in argon 
and krypton, according to which the widths for the counters used in the present 
€xperiment would be expected to be 0:53 and 0-7 respectively. * 

These large fluctuations would present no real difficulty if an unlimited 
number of observations were available. When, however, there are only 
comparatively few observations, ordinary statistical fluctuations are superimposed 
on those resulting from the physical nature of the ionization process. It is then 
necessary to adopt a fairly elaborate method of treating the results. 

The results of the present experiment were first grouped into seven 
momentum intervals, and a statistical procedure (summarized below) was then 
used for fitting curves to the experimental ionization distributions. This 
enabled us to make use of all the observations, duly weighted, including those 
in the ‘high-energy tail’ at some distance from the peak, and to quote values 
for the most probable ionization and the r.m.s. error which were statistically 
valid. ‘The curve which was fitted to the observations was based on Rothwell’s 
curve for krypton, modified by varying the width (without otherwise changing 
the shape) until the best fit was obtained. With the statistical procedure which 
was adopted the dependence of the computed value of the most probable 
ionization on the width of the fitted curve was very slight. For instance, 
changing the width by 20% altered the computed value by only 2%. The width 
found for the distribution was, in any case, not significantly different from that 
quoted for electrons in krypton. 

It would have been possible, of course, to have chosen any one of a number 
of analytical formulae representing single-peaked skew curves with a ‘ high-energy 
tail’, but the curve for electrons seemed physically a more reasonable choice. 

In treating the results no use was made of the fact that two ionization 
measurements were made simultaneously on each particle. ‘The combination 
of these pairs of observations into some kind of mean is useful only when more 
than one kind of particle is being observed. When, as in the present case, the 
particles are very nearly all of one kind, such a procedure does not give any 
increased accuracy. It has therefore been assumed that the observations can 
be treated as though they were independent measurements on two identical 
particles. For this particular experiment the only advantage of the double 
counter was that information was collected twice as quickly. 


4.4. Summary of the Method of Curve-fitting 

The following is a summary of Behrens’ (1951) method of curve-fitting. 
We start with a large number 7 of observed values of a function A (in this case 
the energy-loss), together with a probability curve ¢, as a function of 
x=(A—A,)/€, which it is believed should correspond to the probability 
distribution of the A’s. ‘The purpose of the method is to calculate the best-fitting 
values of the parameters A, and &, and also the r.m.s. errors in these values. 

In the present case A, is equal to the most probable value of the ionization 
and € determines the width of the distribution. 

Let-(Ay 5, ey Ay. «> Aney, 4,) xepresent the set of observed values, 
and define 4; =¢[(A; — A )/E]. 


* It therefore now seems that krypton was an unsuitable choice for this experiment. This is 
confirmed by Cranshaw (1951), who used an argon counter in an essentially similar experiment, and 
obtained a narrower distribution. 
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If we fix A, and €, the probability of any given value of A arising will be _ 


proportional to [¢(x)|/é. The probability P of any particular combination of 


A’s arising is therefore 


Por II (43) e" 
Then log P=constant—m log €+ 2 log 4;; 
hence 56 = “ |e x EP GE a Site ae ied (1) 
where ¢,’ denotes dé/dx at x =«;. 
Similarly, a aS =f (UES oat ae Tienes (2) 


The most probable values of € and A, are found by making both Y and Z 
simultaneously equal to zero. These two auxiliary functions, defined in 
eqns. (1) and (2), are first tabulated for various values of the two parameters. 
The best values of € and Ay are then found by interpolation. This stage of the 
calculation is very simple in practice, since Y depends sensitively on €, and 
hardly at all on Ao, while the reverse is true for Z. The mean square error in Ay 
can be shown to be £,?=(0Z/0A,) +. The differentiation must be performed 
at a point in the (€, Ay) plane corresponding to Y=Z=0, and along a direction 
defined by Y=0. This expression neglects any contribution from third and 
higher order terms. 

Similarly, the mean square error in € is given by (OY/0€)-!, evaluated at the 
same point along the contour in the (€, Aj) plane defined by Z=0. 

The above treatment is correct only when there are negligible observational 
errors. ‘I'he experimental errors can be taken into account by writing the total 
uncertainty E in the form E=[EF,?+(RE,?/n)}/?, where nm is the number of 
observations, E, is the uncertainty arising from the instrumental error in an 
individual observation, and k is a constant, greater than unity, whose value 
depends on the shape of the distribution 4(x) and on how the errors change 
with A. (For a constant error and gaussian distribution k is equal to one; for 
the conditions of this experiment it is slightly greater.) 

The quantity EL has been evaluated numerically for three of the seven 
momentum groups. In no case does the inclusion of the second term introduce 
an additional error of more than 0-2°%, which is quite negligible in comparison 
with the error in curve-fitting. It is obvious by inspection that similar results 
will hold for the other groups. 


§5. STATEMENT OF EXPERIMENTAL RESULTS 
5.1. Most Probable Ionization of Cosmic-Ray Mesons as a Function of Momentum 


The most probable value of the ionization has been computed for each of the 
seven momentum groups: Group I (1-7 to 2:5 x 10% ev/c) contains mesons of 
momentum less than that corresponding to the theoretically predicted minimum 
of ionization (for convenience we shall refer to such mesons as being ‘below’ 
the minimum). In this region the ionization varies rapidly—it is very roughly 
proportional to the inverse of the momentum—so that the group has been 
made as narrow as the rather small number of observations permitted. 
Group II (2:5 to 3:5 x 10%ev/c), and Group III (3-5 to 5-0 x 108ev/c), both 


Ionization by Energetic Cosmic-Ray Particles 445 


contain approximately minimally ionizing mesons. Groups IV, V and VI 
(5:0 to 9-9, 9-9 to 19-5 and 19-5 to 50 x 108 ev/c respectively), being above the 
proton cut-off, are restricted to negative particles. ‘The momentum intervals 
were chosen to give approximately the same statistical weight to each group. 
Group VII contains all the particles of momentum greater than 5 x 10% ev/c. 
Figure 3 shows the curves which were fitted to the experimental data for each 
of the seven groups. The data are summarized in the accompanying histograms. 


a 
20 Ay Group I +ve and -ve a Group V me 
; eS 0-17 to 0-25 x102 ev/c 0:99 to I-95 * 10" ev/c 
: 0 
20 oe Group II +Vve = -ve 
t / NS 0:25 to 0-35 x10 ev/c i ean ul re 


1:95 to 5 x 10° ev/c 
ie | 


Group I +ve and -ve 
0:35 to 0:5 x10°ev/c 


20 


Number of Particles in each Histogram Rectangle 


0 Pye 
Jet 30 Group VII +ve and -ve 
20 Group ¥ src >5 x 10° ev/c 
0:5 to 0-99 x 10% ev/c 
0 | 2 3 4 5 6 i 0 | 2 3 4 5 6 if 
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Fig. 3. Variation of pulse-height distribution with momentum. 


Momentum Intervals and Sign of Particles 


Ag (kev cm*) 


Momentum (10° ev/c) 


Fig. 4. Vaiues of most probable ionization A, derived from statistical treatment of results. 
p, #, 7: Curves calculated from Landau’s theory for protons, and pu- and z-mesons. 
%-: The arrow on the last point indicates that it refers to all particles of momentum 
greater than 5 x 10° ev/c. 


It is emphasized that these histograms have an illustrative value only; they 
were not used in the fitting of the curves to the observations, which was done 
directly. The histograms are included only because the large diagram, in which 
the observed values of energy-loss were plotted against the observed values of 
momentum, is not suitable for reproduction on a reduced scale. ‘The smail 
arrows in fig. 3 show the most probable ionization for y-mesons, and in one case, 
for protons, calculated from Landau’s formula. 
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In fig. 4 the computed values of the most probable ionization are plotted. 


points. For all except the highest momentum group the mean momentum has 


against momentum; r.m.s. errors are indicated by the vertical extension of the © 


>) 


been chosen as abscissa; the unresolved group, however, has been plotted at : 


5 x 10®ev/c with an arrow indicating that it refers to all particles of momentum 
greater than this value. The curves included in this figure show the most 
probable ionization, according to Landau, for jz-mesons, 7-mesons and protons. 


5.2. The Width of the Distribution for Minimally Ionizing Mesons 
Besides giving the most probable ionization, the statistical procedure enabled 
the most probable width of the distribution to be determined. After taking the 
average of the results for Groups II and III (both containing mesons which 
should be approximately minimally ionizing), and correcting for errors, we find 
Full width of distribution at half-height 
Most probable value of specific ionization 


= 75 +04: 


This value agrees, within the experimental error, with the results found by 
Rothwell for 2 Mev electrons and similar counter conditions. 


5.3. Information obtained from Counter Tray D 

The events which have been described as being of the type (A+ B+ C—D), 
and which therefore appear to be associated with the stopping of particles in 
the bottom 8cm of lead, fall into two categories. ‘The first comprises particles. 
of momenta less than about 2-5 x 10% ev/c, the momentum below which no 
y-mesons should be able to traverse the absorber. For this group the majority 
(80°) of all the particles which would have been expected to stop in the absorber 
were recorded as actually having done so. ‘The agreement is not quantitative 
-because, in the restricted space between the supports of the electromagnet coils, 
it was possible to install only a rather inadequate counter assembly. 

The second category contains negative particles with momenta between 
2-5 =x 108 and 5x10%ev/c, and positive particles in the band 2:5 x108 to 
5 x 10% ev/c. The particles recorded as ‘stopped’ are distributed more or less 
at random over the whole momentum spectrum, and are interpreted as particles 
which have suffered catastrophic absorption in the 8cm of lead, without having 
given rise to penetrating secondaries. ‘hey comprise only 5°, of all particles 
observed within the same momentum limits. 

This information is used in the following section to estimate the fraction of 
m-mesons present in the cosmic-ray beam. 


§6. DISCUSSION OF RESULTS 
6.1. The Proportion of Particles other than u-Mesons contained in Groups I-V 1 

The statement has been made earlier in this paper that virtually all the 
particles observed were p-mesons. The arguments that lead to this conclusion 
are as follows. Protons and other heavy positive particles can be excluded 
because of the method of selection (§4.1). None of the recently discovered 
heavy mesons is sufficiently numerous to be worth considering, and no electrons. 
can penetrate the 15cm of lead placed above the apparatus, so that it is only 
necessary to show that 7-mesons were present in small numbers. 

Fast 7-mesons have an interaction length in lead which is about 160 g/cm?. 
Further, for momenta less than 10°ev/c, they give rise to stars with secondaries 
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capable of penetrating a few centimetres of lead in only about 20% of all 
interactions (Camerini et al. 1951). Combining these two facts, we see that 
35% of all incident 7-mesons (at least below 10° ev/c) should be counted as 
(A+B+C-—D) events. From the figures quoted in §5.3 it can be concluded 
that above 2-5 x 10% ev/c a-mesons do not comprise more than about 13% of 
the total meson component. A very similar conclusion can be reached from 
rough calculations based on the known properties of stars, if it is assumed that 
only the stars produced in the 15cm lead filter contribute to the 7-meson 
component. 


6.2. Results obtained for the Unresolved Group (p>5 x 10° ev/c) 

Much of the value of the results would be lost if the large value of the most 
probable ionization found for the unresolved group could be given some 
explanation unconnected with the logarithmic increase, such as, for instance, 
the presence of heavy particles below the ionization minimum. Quite apart 
from the fact that such heavy particles are known to be rare in comparison with 
mesons, it is fortunately very easy to show that there are no particles which have 
the properties required for such an explanation. Alpha-particles, having a 
double charge, and therefore a four times greater ionization at the minimum, 
cannot contribute to the observations which lie in the peak of the ionization 
distribution. The ionization properties of protons, deuterons and tritons, 
below the minimum, are sufficiently well established for it to be certain that, at 
momenta greater than 5 x 10% ev/c, they would have most probable ionizations 
below the value observed for Group VII. Including them in this group could 
only depress, and not raise, the observed value. We are therefore forced to the 
conclusion that for some, or all, of the particles comprising this group the 
logarithmic increase exists. 


6.3. Statistical Significance of the Results: Conclusion 

To estimate the significance of the results we have used a slightly modified 
form of the yx? test. This test deals with any observations in which the numerical 
value of the r.m.s. error is equal to the square root of the quantity observed. 
The tables can therefore be used directly if the results are first scaled so that this 
condition is satisfied. ‘The test has been used in this manner to check the 
agreement of the values of the most probable ionization, computed for 
Groups II to VI, with the predicted curves for p- and 7-mesons. (It did not seem 
desirable to include Group VII because it comprised a mixture of particles of 
unknown momentum and doubtful identity.) 

It is usually assumed that satisfactory agreement between experiment and 
hypothesis has been shown when the probability values derived from the test 
lie within the limits 0-3 and 0-8. ‘The x? value obtained for the u-meson curve 
was 4:3, and for the z-meson curve 2-6, corresponding (for five degrees of 
freedom) to probabilities of 0-5 and 0-75 respectively. ‘The experimental values 
can therefore be satisfactorily accounted for on the basis of Landau’s theory 
although the results are insufficiently accurate to distinguish between the two 
kinds of meson. ‘The evidence given earlier that the particles consisted 
predominantly of p-mesons is nevertheless still relevant, since it has been 
suggested that the absence of a logarithmic rise might be a specific property of 
these particles. 
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A test has also been made to see whether the results could equally well be in 
agreement with the hypothesis that no logarithmic increase exists. No constant 
value of ionization could be found which had a chance of more than one in 10° 
of being in agreement with the results from Groups II-VII. 

The results therefore suggest very strongly that a rise similar to the predicted 
logarithmic increase exists, and therefore disagree with the findings of 
Goodman et al. (1951). There is no obvious reason why there should be this 
contradiction between two very similar experiments, but in view of the 
disagreement we are repeating the measurements with substantially increased 
accuracy. 
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LETTERS TO THE EDITOR 


**Kr—The Possible Existence of a Nuclear Level preceding 
the Metastable State 


During cloud-chamber studies of the angular distribution of successive conversion 
electrons in radioactive isotopes of bromine, many photographs have been obtained showing 
a fast and a slow electron with a common origin in circumstances where the only bromine 
isotope present was **Br. This isotope is known to undergo f-decay to **Kr, followed 
by a 32-7 kev gamma-ray from an excited state of half-life 1-9 hours; another gamma-ray 
of 46 kev with a similar half-life has been reported by Helmholz (1941) but has not been 
found by Bergstrém, Thulin and Andersson (1950) who made spectrometer measurements 
with **Kr separated from fission products. 

The tracks of the fast electrons in the magnetic field applied to the chamber show a 
spread of initial curvatures too wide to be explained as scattering of high-energy 
conversion-electrons, and they must therefore be due to f-particles. Since the resolving 
time of the chamber is about = sec, one particle must follow the other within this time 


100 cs 
if they are to be observed as a pair of tracks of similar width. No f-transition of about 


83 Br 83 Kp 


yo or 2 - 


da = < 10sec A 


46 kev E.2 or M.i 


Possible decay scheme for **Kr. 


1 mev energy (known for **Br) can have so short a half-life, and hence the low-energy 
member of the pair must be the later one and is most probably the conversion electron 
produced in the decay of a short-lived excited state in *“*Kr. Liberation of an atomic 
electron by the f-particle (Bruner 1951) seems too unlikely to explain the number of pairs 
which form a substantial part of all fast electron events observed with “Br. Of the possible 
positions of the short-lived level, that shown in the decay scheme at A (see figure) seems 
most consistent with the allowed f-transition, and the proposed spins and parities of the 
other states as taken from Goldhaber and Sunyar (1951) and from Feenberg and Hammack 
(1949). 

The y-ray from state A has been given the definite energy of 46 kev because of Helmholz’s 
spectrometer measurements, and it is consistent with rough measurements of curvature 
which have been made on the slow electrons in the cloud-chamber photographs. This 
decay scheme seems to explain the results of Helmholz and of Bergstrém et al. The 
failure of the latter to detect the 46 kev gamma-ray would be expected owing to the rapid 
decay of state A during the transference of the source from the separator to the spectro- 
meter. The long half-life reported by Helmholz could be due to the parent B-decay from 
838Br (2-2 hours) because the source used was a film of selenium which had been bombarded 
with deuterons and then placed directly into the spectrometer. In this source **Br can 
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be produced directly or indirectly as shown, and the indirect production by the decay 
of ®8Se could explain also the initial constancy of activity observed by Helmholz. 


82Se (d, n) —————>  *Br ————>"Kr 
22h 
t/~ 
Jo Sis 
VY 


weed p) Seaeaes 


In the present experiments deuteron-bombardment of selenium was again used to 
produce **Br, but the bromine was separated chemically before insertion into the cloud- 
chamber in the form of ethyl bromide. 

On the scheme given the 46 kev transition could be of either electric quadrupole or 
magnetic dipole type, for which respectively a half-life of 10~-* sec or 10~* sec is expected 
if the transition probabilities of Weisskopf (1951) are corrected for internal conversion 
by means of the data of Hebb and Nelson (1940). The longer half-life makes the electric 
transition particularly interesting and it is compatible with the K/L conversion ratio 
reported by Helmholz, whereas the magnetic transition is not. In view of this it seems 
valuable to get more accurate information on the K/L ratio, and experiments have been 
started in conjunction with Dr. P. T. Barrett, who is also searching for short half-lives in 
the krypton transitions. 

I am grateful to Dr. W. I. B. Smith and the cyclotron team for their help with the 
bombardments, and to Professors P. B. Moon and W. E. Burcham for valuable discussions. 
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University of Birmingham. 
1st April 1952. 
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**Kr—An Investigation of the Radiations succeeding 
the Metastable State 


Measurements by Bergstrém (1951) with a beta-ray spectrometer upon electro- 
magnetically separated *°Kr™ have shown that the radiations from the 114 min metastable 
state consist of two highly converted gamma-rays of energies 32:2 kev and 9:3 key. 
Proportional counter examination of these radiations shows that the transitions are in 
cascade and that the lifetime of the 9 kev level is less than 10~’ sec. 

The **K1™ was formed as the decay product of *°Se and **Br produced by bombardment 
of metallic selenium with 7 Mev deuterons; the krypton was separated by melting the 
selenium in a side-tube to the counter, whilst free bromine was retainedin a liquid-air 
trap using SeBr, to provide the carrier. 

The pulse-size distribution was recorded on moving film and analysed by projection in 
a micro-film reader. The peaks of the spectrum occurred in the following manner: at 
9 kev when the 32 kev gamma-ray escaped detection, at 41 kev when both transitions were 
detected in coincidence (9+32=41), and at 28 kev when K-conversion of the 32 key 
transition was followed by a krypton K x-ray of 13 kev energy having an escape 
probability of approximately 99°94, (9 4+-32—13=28). 

From the asymmetry of the peak at 28 kev it is possible to infer the existence of an 
unresolved 32 kev line which would occur when the 9 kev gamma-ray escaped from the 
counter. It is conceivable that the small: prominence at 19 kev was, in fact, the 
corresponding escape peak (32—13=19). 
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‘Two additional pulse-size distributions of the higher energy peaks recorded with 
resolving times of approximately 10-4 and 10~’ sec showed little or no difference when 
normalized at an energy of 41:5 kev. On this evidence it seems probable that the half-life 
of the intermediate state is less than 10~’ sec. 

The ratio of the K to L conversion coefficients of an isomeric transition can be 
calculated from a proportional counter spectrum, providing the fluorescent yield of the 
isotope and the escape probability of the K and L x-rays are known. In this instance, the 
L and M x-rays of krypton were detected with 100% efficiency, whilst the escape probability 
of the K x-rays was very nearly 100%. 'The fluorescent yield for krypton is 58°% for the 
K shell; hence, the total area under the peaks at 19 and 28 kev corresponds to 0:58 Nx, 
where Nx represents the number of K conversion electrons emitted in the 32 kev 
transition. Similarly, the 32 and 41 kev peaks correspond to 0:42 Nx+Nyim. This 
reasoning leads to a value for the K/L+M conversion ratio of 0:32, which is in good 
agreement with that found by Bergstrém. 

Goldhaber and Sunyar (1951) propose that the 32 kev transition of ®*Kr is an electric 
octopole leading to a low lying state of spin 7/2 and even parity. The 9 kev transition to 
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the ground state with spin 9/2 involves a spin change of unity; if the parity changes, it is 
likely to be an electric dipole; if not, a magnetic dipole with a possible admixture of 
electric quadrupole. The values for the corresponding total internal conversion coefficients 
shown in the table were derived from the non-relativistic formulae of Hebb and Nelson 
(1940). These figures were used to correct for internal conversion in Weisskopf’s 
relationship for the expected half-life T. 


Elec. Dipole Mag. Dipole Elec. Quadrupole 
N,/N, Tz 8 PES AN 
T (sec) 10-11 5x 10-8 10-5 


Since the 32 kev transition in **Kr is not totally internally converted, it is impossible 
to calculate an accurate value for the conversion coefficient of the 9 kev transition from 
the proportional counter spectrum; an estimate of this may be obtained as follows: 

N,N. No. of converted 9 kev transitions 
e'""”~ No. of unconverted 9 kev transitions * 
Ty + Log+ Lar 
Thy + Lg2 
where J corresponds to the area under the ‘ line’. 

On the calculated conversion coefficient of 10+5, and until a more accurate 

measurement of N,/N, can be made, the nature of the 9 kev transition can be tentatively 
30-2 


Therefore, NaN ye =10+5 
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identified as magnetic dipole; the isomeric levels of **Kr would therefore have the spins 
and parities shown by Walker (1952). 

The author is grateful to the Department of Scientific and Industrial Research for a 
maintenance grant and to Professor P. B. Moon for many stimulating discussions. 


Department of Physics, P. T. BARRETT. 
University of Birmingham. 
1st April 1952. 
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A Note on Time Reversal in Polarized Nuclear Processes 


It has been demonstrated by various authors (e.g. Wolfenstein 1949, Blin-Stoyle 1951) 
that in a nuclear process in which the bombarding particles (spin 4) are polarized, the 
angular distribution of the outgoing particles in the centre of mass coordinate system has 
the form 

SO)=1O[1 +n. Pf(Op— ws ee (1) 
where n is a unit vector normal to the plane of the process and J(6) is the angular distribution 
obtained when the polarization P; is zero. 

Conversely, if a nuclear process in which the initial particles are unpolarized results 
in outgoing polarized particles of spin 3, their polarization Py has the form P)=ng(6@). 

In a recent paper (Dalitz 1952) the particular case of elastic nucleon-nucleon scattering 
is considered and it is shown that provided the Hamiltonian is invariant under rotation 
and reflection of the coordinate system and time reversal, then f(@)=g(6). 

It is the object of the present note to show that this result is a particular case of a more 
general theorem. We shall, in fact, prove for the inverse nuclear processes 


A+P—+B-+Q (denoted by I) and B+Q--A-+P (denoted by II) 
in which the particles P have spin } and are polarized, that 


FiO)=2n@).. "eee (2) 


(Elastic scattering is implied, of course, when A and P are identical with B and Q 
respectively.) 

Referring to a previous paper (Blin-Stoyle 1951) it can be seen by slight rearrangement 
of the results given there, that f{(@) and gi;(@) are given by the expressions 


= |2(—1)4875"0) Yi Ol*f(6)=2i » {a(— 14 Spb“ (D YL (@)}* 
aog 


abpq Ll 
x {a(— 1)ES a (DV 9000) | nor eee eee (3) 
and 
= |5(—1)Sehe (IL) V7 (4) |2er1(8)=2i & {E (— 1452 (I) Y¥()}* 
abpq Li abq Ll 
x {2( ={is aie 6), ee (4) 
where the scattering matrix Ry enl 8) (cf. Wigner and Eisenbud 1947) is related to the 
function Fone (1) defined in the above paper by 
.apLM i; i I 
Shain (1) =t-L21- Pi )MeR P(g ae a pee ee (5) 


with an analogous relation for Seer (il). 


a, p, b and q refer to the z-components of the spins of A, P, B and Q, and L. M and 
/, m are the orbital angular momentum quantum numbers representing the relative motions 
of A, P and B, Q respectively. The factor Sopa Meee implies the conservation of 


the z-component of angular momentum and removes the necessity of summing over 
M and M’ in (3) and (4). . 
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It is also demonstrated in an appendix to the above paper that 
LM arp i= 
Semen Eo hime kimmel ae (6) 
only one of the signs holding for a given process. 


Now, if the Hamiltonian is invariant under reflection of coordinate axes and time 
reversal, it can be shown (e.g. Coester 1951) that 


Sapiall)= + S—3—gi me (1)= + Som” (1) 
by (5) and (6), the relevant sign being determined by the process under consideration. 
It therefore follows at once from (3) and (4) and the fact that L+/ is always even or 
always odd for a given process, that the equality (2) is true provided the Hamiltonian 
possesses the invariance properties listed. 
I am indebted to Professor Pryce for informing me that the foregoing theorem is an 


example of a more general relation which can be shown to hold between inverse physical 
processes. 


The Clarendon Laboratory, R. J. BLIN-STOYLE. 
Oxford. 
19th March 1952. 
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Neutron Flux and Spectrum Measurement 


The measurement of a neutron spectrum by observing the recoil protons produced in 
a hydrogenous medium is usually subject to two difficulties. If the energy spectrum of 
all the protons produced in an ion chamber containing hydrogen is to be measured, the 
dimensions of the chamber must be large compared to the range of a recoil proton. For 
protons of more than a few hundred kilovolts energy this requires either a very high 
pressure filling or a very large chamber. If, alternatively, the proton recoils from a thin 
hydrogen-containing film are observed, then, although the chamber dimensions may be 
reduced to one proton range, the number of recoils from a thin film is inconveniently small. 
An attempt to overcome the first difficulty is described here, in which a scintillating 
anthracene crystal is used instead of an ion chamber. A crystal of quite small size 
(~5 mm side) will fulfil the condition as to proton range, and, provided the relation 
between light output and proton energy is known, it should be possible to determine the 
energy spectrum of proton recoils in such a crystal. Scintillations from the recoiling 
carbon atoms are negligibly small since not only is their energy lower, but also the light 
output per Mev dissipated in the crystal is much less. The relation between light output 
and proton energy has been calculated by Birks (1951) in terms of the energy—range relation 
for protons, and is given as 
dS A 
dE 1+BdE/dx’ 


where S represents light output, A and B are constants, and dE/dx is the rate of energy 
loss of a proton in anthracene. 

Anthracene was chosen for this work because of its high luminous efficiency, which is 
important since the statistical spread in the number of electrons released at the 
photocathode must be kept as low as possible to give the best energy resolution. 
Unfortunately, the use of an organic scintillator gives the counter a disproportionately 
high sensitivity to y-radiation, and its use is limited to measurements in the absence of 
y-radiation, or where the y-radiation effect may be subtracted by means of a ‘background’ 
run. 

A counter consisting of a clear piece of anthracene, weighing 0-433 g, cemented to the 
window of an E.M.I. photomultiplier type VX 5045 was exposed to 2:54 Mev neutrons 


Counts per Channel 


454 Letters to the Editor 
emitted from the D-D reaction at 90° to the beam, and a pulse height distribution was 
taken with a 30-channel pulse analyser (fig. 1). Using eqn. (1), this distribution was 
transformed to the energy spectrum of the recoil protons to give fig. 2. On the same 
figure is plotted the neutron spectrum N(E) derived from 


1H) ye) 
o(E) dE 
N’(E) is the distribution of proton energies, and y(E) is the n—p scattering cross section. 


To test the counter further it was used for an absolute measurement of the flux from 
the D-D source, under conditions for which this flux was known. ‘To this end the counter 
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and an Ilford type C2 nuclear plate were simultaneously exposed at equal distances from 
the source, and the number of pulses above 27:5 volts amplitude was counted. From this, 
and from eqn. (1), the total number of proton recoils was calculated, and so the number 
of neutrons falling on the crystal. This figure was compared with that obtained from 
examination of the plate, and the results are: crystal 8:7(+0-2) x10’ neutrons/cm?, 
plate 8:0(+0-5) x10’ neutrons/cm?. 


Atomic Energy Research Establishment, M. J. Poo.e. 
Harwell, Berks. 
7th March 1952. 
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The Range-Energy Relation for Protons in Aluminium 


No experimental determinations of the stopping power of aluminium for protons in the 
energy range from 5 to 10 Mev appear to have been reported. ‘To fill this gap, Ilford C2 
nuclear emulsions have been used to measure the energy loss of initially 10 Mev protons in 
traversing a thickness of aluminium which reduces their energy to about 3 Mev. 

The source of protons was the extracted beam of the Birmingham 60 in. cyclotron, 
operated with no ion-source arc to.give a suitably low intensity. The plates were tilted 
slightly from the horizontal so that the protons entered the emulsion at a small glancing angle. 
Part of the beam struck the plate after passing through a sheet of aluminium mounted just 
before the plate with its plane vertical: part of the beam passed below the aluminium and 
struck the plate directly. The difference in height of the two beams is only about a milli- 
metre, anid is in the vertical direction in which the beam energy is known not to vary greatly. 

The table shows the results of six determinations, including measurements made with 
four different plates and three different absorbers. 20 to 30 tracks were measured at each 
range, the mean range being thus determined within 2. In this table, U refers to the 
unstopped beam, S to the beam which has traversed the aluminium absorber. Columns 4 
and 5 are obtained by converting ranges in the emulsion into proton energies by means of the 
range—energy relation for C2 emulsions given by Rotblat (1951), and then converting these 
energies into ranges in aluminium by means of the range-energy relation calculated by Smith 
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(1947). Column 6, which is the difference between columns 4 and 5, is the thickness of 
aluminium which the beam would have traversed if Smith’s range-energy relation had been 
accurate. Column 7 gives the actual thickness of aluminium traversed. : 

Column 8 lists, as a percentage of the absorber thickness, the difference between the 
known thickness of aluminium and that predicted from Smith’s range-energy relation. The 
errors quoted are probable errors calculated from the spread of the measured track lengths 
in the emulsion, and do not include the uncertainty in the range-energy relation for the 
emulsion. 

In all cases, Smith’s figures predict too low a value for the thickness of aluminium. In 
column 8, the agreement between the average for plate 1, (1:0+0-2)%, and the average for 
the other three plates, (0-6+0-2)%, suggests that there is no serious systematic difference 
between plates. The final average of column 8 suggests that Smith’s range figures are too 
low by (0:8+0-14)% ; when the $°% uncertainty in Rotblat’s range-energy relation is 
included, the error in Smith’s ranges becomes (0:8+0-5)%. 'This correction applies 
essentially to the range at 10 Mev, since the range of the stopped beam is so much smaller 
than that of the unstopped beam. 

The stopping power of a substance for heavy particles may be written (Smith 1947) 


452 an 
EON Ie a tenis eee (1) 


where — dE/dx is the rate of energy loss of the particle, ez the charge of the particle, v the 
velocity of the particle, Z the atomic number of the stopping substance, N the number of 


1 ? 3 4 5 6 7 8 
4 : : 2, Thickness of Al Mera 
Plate Range in Emulsion (j1) Range in Al (mg/cm?) traversed (mg/cm?) Mensured sath (%) 
U S U S Smith Measured MSE 
1 498-9 41-6 sli} 11-4 (39:8 mel 412 1-:0+0-4 
i 492-6 38-2 149-4 10:5 138-9 141-2 1-640-4 
1 490-8 30-0 148-6 8-0 140-6 141-2 0-4+0-2 
2 512-9 55-0 155-6 15e2 140-4 = 140-6 0-2+0-4 
3 513-2 91-4 155-7 25°8 129-9 131-3 1-1-0-4 
4 524-2 70-1 159-4 19-6 139-8 140-6 0-6-+0:2 


stopping atoms per unit volume, J the average excitation potential of the electrons of the 
stopping substance; e, m are the charge and mass of the electron. Ranges are obtained 
from (1) by numerical integration of the expression 


EB dE 
R= | ae ene (2) 


The parameter J cannot be calculated accurately and the value of 150 ev used by Smith is an 
empirical one (Wilson 1941) on which the probable error is more than large enough to account 
for a 1% range error at 10 Mev. If the error in Smith’s ranges is assumed to be due entirely 
to an error in J, then it may be shown that the percentage range error should vary little over 
the range of proton energies from 10 to 80 Mey. 

Bloembergen and van Heerden (1951) have measured the ranges of fast protons in 
aluminium, and have found, also, that Smith’s ranges are too low. ‘The corrections they 
find are (0-9 + 0-4)% between 56 and 76 mev, (1-4 +0-9)% between 35 and 52 Mev. These 
results, and that quoted above for 10 Mev protons, are consistent with the assumption that 
the error.in Smith’s ranges is due entirely to an error in his value of J. 

From eqns. (1) and (2) it may be shown that 
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where AR is the error in range caused by an error AJin J. ‘The present experiment implies 
a value J=155 + 3 ev, and is to be compared with 159+ 5 ev found by Bloembergen and van 
Heerden at energies between 35 and 76 Mev, and with 156+ 3 ev deduced from Sachs and 
Richardson’s (1951) direct measurement of the stopping power of aluminium for 18 Mev 
protons. The three experimental values are in good agreement, and yield a mean value 
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I=156+2ev. This figure has been used in conjunction with eqn. (3) to prepare the curve, 
which gives the percentage correction which should be added to Smith’s tabulated ranges 
from 2 to 80 Mev. ‘The range—energy relation corrected in this way is unlikely to be in error 
by more than 4%. 
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Percentage correction to be added to the proton ranges in aluminium tabulated by Smith (1947). 


The author’s thanks are gratefully offered to Dr. J. Fremlin for helpful discussion, to 
Mrs. E. Munday for her aid in processing plates and measuring track lengths, to Mr. W. 
Hardy who operated the cyclotron, and to the University of South Africa for the award of a 
scholarship. 


Department of Physics, D. H. Stumons. 
University of Birmingham. 
18th January 1952. 
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The Response of a Nal Scintillation Counter 
to Slow and Fast Neutrons 


The fluorescence of NalI(Tl) crystals under charged particle and y-ray excitation has 
been described by Hofstadter (1950) and others. We have used a neutron source, free 
from y-ray background, to determine its response to slow and fast neutrons with a view to 
using it to investigate the inelastic scattering of fast neutrons. 

An E.M.I. photomultiplier fitted with a crystal of Nal(Tl) (143x108 mm) was 
mounted close to the target of a D-+D neutron source operated with a 40 wa resolved beam 
of 190 kev deuterons; the target was formed by the occlusion of deuterium in an 
aluminium sheet 0-3 mm thick which gives rise to a negligible y-ray background (Beghian 
et al. 1950). Slow neutrons were produced in a tank of water surrounding the target, 
whereas for measurements with fast neutrons the crystal was placed in contact with the 
target. 

The y-rays produced by the capture of slow neutrons in iodine gave a continuous 
pulse height distribution extending up to about 6:2 Mev. This limit probably corresponds 
to a y-ray of 7:2 Mev detected by pair production in the crystal (cf. 7-0 Mev, Kubitschek 
and Dancoff 1949). No distinctive peaks can be identified in the spectrum although a 
steep rise in the vicinity of 2-2 Mev can probably be attributed to neutron capture radiation 
by the hydrogenous moderator. Intense 478 kev y-radiation arising from 1B(n, «)?Li* 
in the glass of the phototube was minimized by mounting the scintillation crystal on a 
‘ Perspex ’ pillar. The efficiency of the crystal for the detection of y-rays arising within 
it was estimated from the photoelectric, Compton and pair production cross sections, and 
this was used to assess the number of quanta originating in the capture process; the number 
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of capture processes was determined from the B radioactivity (1°81) in the crystal. In this 
way the number of y-rays with energy greater than 3 Mev was found to be of the order of 
1 quantum per neutron captured. 

The pulse height distribution due to fast neutrons cuts off at about 2-5 Mev. We 
attribute it to y-rays due to inelastic scattering. At 450 kev there is indication of a complex 
photoelectric line which must be attributed to a number of y-rays of closely similar 
energy : although recoiling sodium nuclei resulting from direct collisions with neutrons 
have a continuous spectrum with a maximum energy of 450 kev, it does not seem possible 
that they could account for a peak. 

A chromium disc (2 cm diameter x 0-4 cm thick) was then mounted in contact with 
the target and the pulse height spectrum due to y-rays from inelastic scattering of the 
neutrons in the chromium was obtained : the difference counting rate with and without 
the scatterer is shown in the figure. Although the background due to primary fast 
neutrons is large, there is a good indication of a y-ray at 1-4 Mev in agreement with earlier 
measurements (Grace et al. 1951).. No lower energy y-ray is detectable above 250 kev. 
A rise in the Compton distribution at low energies is seen which is probably due to 
degradation of radiation in the chromium and to the detection of neutrons which have 
suffered inelastic scattering. The geometrical arrangement was such that the presence 
of the chromium scatterer did not interfere with the direct neutron beam from the target 
to the crystal. 

The cross section for this process was estimated to be 1:5 x 10~*4 cm?, in fair agreement 
with an earlier figure (Grace et al. 1951) of 1:2+0-4 x 10-4 cm?. 
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Pulse height distribution due to y-rays from chromium excited by fast neutrons. 


This method of investigating inelastic scattering seems to be limited to nuclei with 
well-spaced levels and large cross sections. It is hoped that with the use of a collimated 
fast neutron beam the method may be made more sensitive. 

We are greatly indebted to Dr. F. V. Price for putting the facilities of the high voltage 
set at our disposal. 


The Clarendon Laboratory, M. A. GRACE. 
Oxford. H. R. Lemmer. 
26th March 1952. H. Hasan. 
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The Kinetic Formulation of Conduction Problems 


Problems in electrical conduction are usually formulated by means of the Boltzmann 
equation. An alternative formulation, using an integrated form of the Boltzmann equation 
derived by simple kinetic theory arguments, has proved useful in solving certain ‘thin- 
conductor’ problems (Chambers 1950, 1951). Recently Huxley (1951) has given a lengthy 
‘free path’ formulation, in which the current density is found by first finding the position of 
the centroid of a group of electrons and then taking its time derivative. His final equations 
are exceedingly complicated. It seems worth while to point out, therefore, that a much 
simpler and more general formulation may be obtained by calculating the current density I 
directly on a kinetic theory basis. 

Consider a volume element dr=dx dy dz at point rp in the conductor, and electrons 
passing through it with energy E and velocities in the range dv=dv, dv, dv, about V9, or with 
wave numbers in the range dk about ky. For free electrons, #k=mv; for electrons in an 
ionic lattice v=(1/h) grad, E. Let fy be the unperturbed electron distribution function 
(so that fy dr dk/477* is the number of electrons in dr, dk), f the perturbed value in the presence 
of applied fields, and 6f=f—f». Let the relaxation time for electron collisions be 7, and 
assume that the distribution function of electrons immediately after collision corresponds to 
6f=0. Now electrons passing through ry with velocity vj) and energy E at time t, will have 
followed a certain trajectory since their last collision; the value of f(r, Vo, to) is found 
simply by integrating the number scattered into the trajectory at previous points along it with 
energy E—AE (where AE is the energy then acquired from the applied fields before 
reaching ro), weighted by their probability of reaching ry : 


ty 
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Hence, on writing f)(E—AF)=f,)(F)— AE df,/dE and integrating by parts, we obtain 


, atl) fl E— AE) expt—(to—2)/7} . 


Bf (to Vor to) = —(dfaldE) |" (di[r) AEexp{—(ty—0/7} 
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to ; : : 
Now AE(t)= —e| E(r,s) . v(s) ds where E (r, t) is the applied electric field (magnetic fields, 
if present, have no direct effect on AE, but alter v(t)) and d AE/dt=eE(r, t) . v(t); further, 
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and lastly, if r=7(v)=7(t), exp{—(t)—1t)/7} must be replaced by exp{— |" ds/t(s)}, so that 
finally we obtain : 


I(r); t))=— | v Teak dt E(r, t) . v(t) ex | ds/t(s)} 1 
a 47° all k dE —o ; ; Ba t aS (1) 

which is the general solution of the Boltzmann equation for any isothermal conductor in 
stationary or varying, uniform or non-uniform fields E and H, provided that the conductor 
is of dimensions large compared with the mean free path /=7|v|=7v, and that the charge 
distribution is uniform, i.e. that [of dk=0. If [df dk40, a simple ‘ diffusion’ term appears; 
and for a general formulation of ‘thin conductor’ problems the equation needs only slight 
modification (cf. Chambers 1950, 1951). More important, since eqn. (1) can be applied to 
electrons in an ionic lattice, where v=(1/f) grad, E, it can be used (Shockley 1950) to 
elucidate the difficult problem of magneto-resistance effects in metals, ; 

For free electrons, assuming t=7(v) only (whether in a metal, semiconductor, or 
ionized gas), eqn. (1) gives, for H=0 and E uniform, 
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Huxley’s equations (1), (38), and (41), for the special case /=rv=constant, yield (2) after 
sufficient integration and rearrangement. For H=H, (uniform and constant) (1) becomes 


(I,+ily), I= 


Sretm2le dd to 
- ait { v° Pa” - dt[(E,(t)+1E,(t)) expiw(to—t), E,(t)] exp{—(to—t)/T} 


0 
Ree (3) 
where w=— eH,/me. ‘This relatively simple expression replaces and generalizes Huxley’s 
equations (1) to (5) and (38). 
Royal Society Mond Laboratory, R. G. CHAMBERS. 


Cambridge. 
20th March 1952. 
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A Combinatorial Problem in Electron—Photon 
Cascade Theory 


In electron—photon cascade theory in approximation A the following problem arises ; 
to find the number of different cascades in which one primary electron or photon gives rise 
to given numbers of electrons and photons, the cascades being distinguished only by 
the order in which the bremsstrahlungs and pair productions occur. A somewhat similar 
combinatorial problem in counting cosmic rays has been considered by Schrédinger (1951). 

Suppose that a primary j gives rise to 2n+j electrons and m photons; for a primary 
electron (j=1) the number of electrons must be odd and for a primary photon (j=2) the 
number even. The numbers of electrons and photons before the last collision are 
2n+j and m—1 if this collision is a bremsstrahlung and 2(n—1)-+-j and m-+-1 if the collision 
is a pair production. Hence, if N(n, m) is the number of possible cascades, 


N(n, m)=N(n, m—1)+N(n—1,m4+1).- - a a ee (1) 
The solution of this difference equation is 
m+1 (2n-+-m)! Q) 
n+m+1n!(n+m)!’ : 
satisfying the obvious initial condition V(0,0)=1. The values of N(n, m) for small values 
of n, m are given in the table. Each entry in the table is found by adding to the number 
on the left the number on the right of the one above. The value of N increases rapidly 


with increasing n and m; the number of different possible cascades in which a primary 
photon gives rise to 20 electrons and 20 photons exceeds 108. 


N(n, m)= 


m\ 0 1 2 3 4 
0 1 1 1 1 1 
1 1 2 3 4 5 
2 2 5 9 14 20 
3 5 14 28 48 75 
4 14 42 90 165 215 


The cascades may be represented by the different ways of ordering m-+-n 1’s and 
n 2’s such that from the left there are never more 2’s than 1’s. A1 stands for a 
bremsstrahlung, a 2 for a pair production, and the orderings represent all collisions if the 
primary is an electron and all except the first (which must be a pair production) if the 
primary is a photon. 

For example, if a primary electron gives rise to five electrons and one photon, 
n=2, m=1, and the different possible cascades (five in number) are represented by 
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11122, 11212, 11221, 12112, 12121. These orderings also represent the cascades in which | 
a primary photon gives rise to six electrons and one photon. The ordering 12121, for 
example, represents the cascades given in the figure. 
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The cascade 12121. 
(a) Primary electron giving five electrons, one photon. 
(6) Primary photon giving six electrons, one photon. 


The author wishes to thank Dr. H. Messel and Mr. G. Szekeres for helpful discussions. 
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The Parameters of Partially Degenerate Semiconductors 


In semiconductor problems, the concept of the Fermi level is of considerable importance, 
and some interest attaches to methods for its evaluation. For this purpose, consider a 
simple n-type semiconductor, with an impurity concentration of Ng donor centres per unit 
volume, each having an energy — ¢g with respect to the bottom of the conduction band. 
Then the free electron concentration ng and the Fermi level e* can be determined from 


Ne=No(Me/ M)?!? Fy) y (n*) Bro eles roll) 
and Ne=Nal(1 + exp (n*+4)]7, eRe) 
where n*=e*/kT, nd=ea/RT, mo=40(2mkT] h?)3/? 

pl?  atdn 
d : a OS Oo 
an Fy) (7 *\= I. ese *) (3) 


Methods of solving these equations have been developed by Shifrin (1944) and Putley (1949). 
amongst others and are discussed in detail elsewhere by the present author (Blakemore 1952). 
The Fermi-Dirac integral F’\,.()*), which enters into the calculation has been accurately 
tabulated for 7*=—4(0-1)20 by McDougall and Stoner (1938). For the evaluation of 
semiconductor parameters, however, it is often useful to have simple analytical expressions 
for this integral. 
In cases of strong degeneracy, a good approximation is 


Fy) 9(q*) 2% *8/? [1 +ar?/8n*?] Ste geste (4) 
which is satisfactory to within a few per cent for 7* greater than 1:5. ‘The range of appli- 
cation can be extended downwards by observing that (4) is identical with the first two terms. 


of the expansion 
Fj 9(n*) = [n*? +-777/6]9/4 ney iE!) 
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which has an error of less than 3°% for n* greater than 0:8, as shown by curve 6 in figs. 1 and 2. 
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Fig. 1. The Fermi-Dirac integral, Fy/2(*), 
[curve a], together with the approxi- 
mations 6 of eqn. (5), and c of eqn. (6). fig. 1. 

For semiconductors which are slightly degenerate, modifications of the classical approxi- 
mation F,,()*) =4371/? exp (y*) which applies for large negative values of »*, have been 
proposed by Ehrenberg (1950), and by Landsberg, Mackay and McRonald (1951). Both 
of the suggested modifications are of the form 

Fy) (*) =4n7!*[6--exp (7 ee 

where Ehrenberg suggests 6=0-25, and Landsberg et al., 6=2-*/?. The latter authors 

point out that different values of ¢ are useful for various small ranges of 7*. By the sub- 

stitution 6=0-27, it can be seen from curve c in figs. 1 and 2 that the function can be 
represented for all negative values of y*, and for positive values up to 1:25, with an error of 


less than 3%, 

The use of Fj2(n*) 47 
F,2(y*) to be covered over the complete range of n*. 
over from one expression to the other is at 7*= +1. 


Fig. 2. Errors introduced by the approxi 
mations given in the curves 6 and c of 


1/2(0-27+exp(—7*)]~! in conjunction with (5), enables 
A convenient poini at which to change 
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Schumann Region Absorption Spectrum of Copper Vapour 


Amongst atomic spectra, that of Cur has occupied a special place in presenting an 
abundance of peculiarities arising from series perturbation and auto-ionization. ‘The 
spectrum has been extensively studied, notably by Shenstone, who has summarized existing 
knowledge of it (1948). Included in Shenstone’s term lists are a number of tentative 
assignments to the configuration 3d°4s5p. These assignments were made from measure- 
ments of anumber of emission lines lying in the 1570-1750 A range, classified as combinations 
with the 3d°4s? 2D metastable term. The only other allowed combinations of the 3d*4s5p 
levels are with the ground state, and are thus expected to lie in the short wavelength part of 


the Schumann region. 
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Partly in view of the tentative nature of Shenstone’s assignments, partly in the hope of 
observing combinations of the ground term with levels of 3d°4s6p, the absorption spectrum 
of copper vapour has been photographed over the range 1060-2500 A. The copper vapour 
was produced in the usual way in a King furnace (Garton 1952), and the hydrogen spectrum 
and the Lyman continuum were used as backgrounds. The plates reveal a set of fifteen 
new lines in the range 1200-1430 A, of which the strongest are a group of six between 1370 
and 1420 A, and six between 1220 and 1300 A. 

Initial studies of these plates indicate that the following of Shenstone’s levels are real 
and combine with the ground state: x*P3)2, x*Pj;2, w’Ps;2. Three of the new absorption 
lines are thus accounted for by upper levels already known. The shorter wavelength lines 
no doubt involve upper levels arising from 3d*4s6p and 3d*4s4f. 

Full details of the spectra, measurements and classifications will be given later. 


Department of Physics, W. R. S. GaRTON. 
Imperial College, London. 
17th April 1952. 
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Further Results inthe General Theory of Barrier Layer Rectifiers, by P.T. LANDSBERG. 


ABSTRACT. It is pointed out that the temperature variation Vy =a—bT of the diffusion 
potential, as deduced from the forward characteristic of rectifiers, is in broad agreement 
with theory. Hence formulae are derived for the temperature variation of the zero resistance 
of a rectifier having an arbitrary distribution of impurity centres, and subject to either diode 
or diffusion theory. Interpretation of experiments with the aid of these formulae leads to 
the temperature dependence of (a) the effective mass of current carriers in germanium 
rectifiers, and (6) of the mobility of current carriers in selenium rectifiers. The latter 
curve has the same form as that obtained from Hall effect measurement, and it is pointed out 
that it is of broadly the same shape as the corresponding curves for high resistivity samples of 
silicon and germanium. It is shown that two d.c. characteristics of the same rectifier, but 
corresponding to different temperatures, can intersect. Formulae for the applied voltage and 
the temperature rise which lead to the onset of thermal instability are deduced for all 
impurity centre distributions, and for diode and diffusion theory. A general theory of current 
creep in rectifiers is outlined, and good agreement obtained with experiments on selenium 
rectifiers. 


Breakdown in Selenium Rectifiers, by R. COOPER. 


ABSTRACT. ‘The paper describes an experimental investigation of breakdown in 
selenium rectifiers. Cooling from about 110°c to —74°c caused the breakdown voltage to 
decrease continuously, and at the same time anomalous behaviour was observed in the 
current—voltage characteristics. When high voltages were applied, the current decreased 
as the temperature increased, but at low voltages rise in temperature caused the current to 
increase. 

The effect of changing the thermal dissipation constant was similar to that caused by a 
change in ambient temperature. Reduction in the dissipation constant, effected by placing 
the rectifier in a vacuum, was observed to increase the breakdown voltage, but the leakage 
current at low voltages also increased. 

None of the available theories appears to account satisfactorily for the observations. 


Characteristics of Scintillation Counters, by G. F. J. GaRtick and G. 'T’. WRIGHT. 


ABSTRACT. Experimental studies have been made of the distribution in size of output 
pulses from a scintillation counter detecting monoenergetic alpha-particles. Separate 
investigations have been made of contributions to the pulse size spread by the photomulti- 
plier and by the scintillating crystal. It is shown that the main effect is due to statistical 
fluctuations in the number of photoelectrons produced at the photomultiplier cathode. 
The effect of the multiplying stages on these fluctuations is investigated and also the effect 
of flaws in the detecting crystal. "The frequency distribution of the pulses is usually found 
to follow the gaussian form. 


The Current in the Electron Immersion Objective, by L. Jacos. 


ABSTRACT. It is shown that for an electron immersion objective under the condition of 
constant cut-off voltage the emission current is uniquely defined by the cross-over potential. 
There is thus some theoretical foundation for the empirically observed law, that the current 
at zero modulation voltage depends only on the three-halves power of the cut-off voltage. 


An Improved Method for Measurements of Optical Constants by Reflection, by 
D> G AVERY, 


ABSTRACT. At wavelengths outside the photographic range the determination of 
the optical constants of absorbing media involves the measurement of intensity changes 
on reflection at the surface of the medium. A review of previous methods for their 
determination is made and the experimental difficulties are indicated. A method is 
described which surmounts some of these difficulties by avoiding the direct measurement 
of reflection coefficients. The accuracy to be expected from the method and an example 
of its use are briefly discussed. 
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Primary Aberrations and Conjugate Change, by C. G. WYNNE. 


ABSTRACT. €xpressions are derived, in terms of the functions normally used in optical 
designing, for the variation with conjugates of the primary aberrations of any symmetrical 
system with isotropic object and image spaces. The possibilities of invariance of the 
aberrations, or their simultaneous correction for more than one magnification, are discussed. 


Propagation of Ultrasonic Waves in Vapours near the Critical Potnt, by 
H. D. Parsroox and E. G. RICHARDSON. 


ABSTRACT. An experimental study is made of the velocity and absorption of ultrasonic 
waves in carbon dioxide and ethylene at frequencies between 4 and 2 Mc/s and at pressures 
up to 100 atmospheres in a variable-path acoustic interferometer. At pressures above 
10 atmospheres the results are independent of frequency (in this range) in the gaseous phase. 
The velocity and absorption coefficient reach a minimum near the condensation point, but 
especially at the critical point itself have high and indeterminate values when the change of 
phase takes place. The ultrasonic viscosities calculated from these measurements are of the 
order of 1000 times those derived from transpiration or low-frequency oscillation experiments. 


The Influence of Electrolytes on the Mechanical Properties of Certain Metal 
Single Crystals, by E. N. pa C. ANDRADE and R. F. Y. RANDALL. 


ABSTRACT. The contact of electrolytes has a marked effect on the mechanical behaviour 
of single crystals of cadmium and zinc which have been exposed to the air. This is due 
to the disruption, or formation, of a thin surface film, ptobably of oxide or hydroxide. 
If ordinary single crystal cadmium wires are loaded so as to flow very slowly, the contact 
of solutions containing free cadmium ions increases the rate of flow : a similar effect exists 
in the case of zinc treated with solutions of zinc salts. If the surface of the wire is cleaned 
by preliminary thermal evaporation there is no such effect, establishing that it is due to 
a surface film normally present. Solutions of cadmium or zinc nitrate stop the flow and 
raise the critical shear stress appreciably, which is attributed to the formation of a surface 
film, probably of hydroxide. With the cadmium nitrate there is a temporary increase 
of flow before the cessation, which is attributed to the cadmium ion, for there is no such 
effect with the zinc nitrate. With polycrystalline wires there is no electrolyte effect. The 
work described supports the view that a thin oxide or hydroxide film increases the 
mechanical resistance of single crystal wires, and that with cadmium and zinc wires such 
a film is normally present. 


PROC. PHYS. SOC. VOL. 65, PT. 6—a (A. LAGERQVIST, G. NILHEDEN AND R. F. BARROW) 
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Fig. 1 (a), (6) and (c). Bands of the p—x system of SiS in absorption. 
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Fig. 1(d). The 2,0 band. Magnification 11 x. 
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